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ABSTRACT 


Network  studies  are  directed  toward  effective  utilization  of  a 
collection  of  seismic  stations  as  a  coherent  worldwide  seismic  network  to 
lower  the  detection  tnreshold  for  seismic  events  and  improve  their  classi¬ 
fication.  This  report  presents  an  analysis  of  ambient  seismic  noise  seen 
at  the  network  level  and  is  directed  toward  characterization  of  th  j  network 
noise  field.  Simultaneously  recorded  noise  samples  from  two  VE^A  stations 
and  eight  LRSM  stations  are  investigated  through  analysis  of  power-density 
spectra,  coherence,  high- re  solution  frequency-wavenumb°r  spectra,  and 
K-line  spectra.  The  ambient  short-period  network  noise  field  is  found  to 
be  highly  variable  with  time  and  with  geographical  position.  Variations  in 
noise  power  are  directly  related  to  regional  storm  activity.  Dominant  noise 
power  at  most  stations  is  trapped-mode  surface -wave  energy  that  does  not 
exhibit  useful  interstation  coherence. 
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CPO 

ACRONYMS 

Cumberland  Plateau  Seismological  Observatory 

LASA 

Large-Aperature  Seismic  Array 

LPC 

Low-pressure  center 

LRSM 

Long-Range  Seismic  Monitor 

MCF 

Multichannel  filter 

SNR 

Signal- to -noise  ratio 

TFO 

Tonto  Forest  Seismological  Observatory 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 

A  partial  goal  of  Contract  F33657-67-C-0708  is  an  examination 
of  methods  to  utilize  effectively  a  collection  of  seismic  stations  as  a  coherent 
seismic  network.  A  subtask  of  the  network  studies  involves  a  comprehensive 
examination  of  broadband  ambient  noise  as  seen  at  the  network  level. 

This  report  presents  the  results  of  the  analysis  of  five  noise 
samples  recorded  at  the  network  stations  between  12  August  1964  and  29 
January  1965.  The  analysis  determines  the  salient  features  of  both  short- 
period  and  long-period  noise  at  the  station  level,  with  particular  attention 
given  to  spectral  analysis,  spatial  organization,  and  identification  of  modal 
content  and  noise  sources.  The  structure  of  the  network  noise  model  is 
then  inferred  from  an  analysis  of  the  combined  effects  of  these  features. 

The  network  consists  of  VELA  arrays  TFO  and  CPO;  LRSM 
arrays  OONW,  GGGR,  NPNT,  HWIS,  and  LZBV;  and  single-instrument  sta¬ 
tions  ADIS,  RKON,  and  DHNY.  This  study,  which  is  restricted  to  vertical- 
component  ambient  noise,  consists  of  single-sensor  analysis  through  power 
spectra  and  modal  noise  analysis  utilizing  array-station  data. 

Outputs  of  single-sensor  stations  (and  single  seismometers 
selected  from  each  array  station)  have  been  studied  to  determine:  the  dis¬ 
tribution  of  station  noise  power  as  a  function  of  frequency,  time,  and  geo¬ 
graphic  position;  the  correlation  between  station  noise  levels  and  regional 
storm  activity;  and  the  seismic  validity  of  the  recorded  data. 

Array-station  data  have  been  studied  to  determine  seismic 
propagation  modes  and  velocities  of  the  coherent  energy,  to  establish  the 
sources  of  this  energy,  and  to  determine,  if  possible,  whether  the  stations 
are  bodywave  noise-limited. 
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Power  spectra  are  computed  from  five  short-neriod  noise  sam¬ 
ples.  Power  spectra  and  interstation  coherences  also  are  computed  from 
long-period  data  for  two  of  the  five  noise  samples. 

Modal  noise  studies  are  conducted  through  the  analysis  of  high - 
resolution  frequency-wavenumber  spectra,  K-line  spectra,  Wiener  disk- 
model  MCF  outputs,  and  straight- sum  outputs. 

Short-period  data  appear  to  be  seismically  valid  and  generally 
usable  between  0.  2  and  2.  0  Hz.  Above  2.  0  Hz,  several  apparently  nonseismic 
noise  peaks  appear.  Generally,  station  array  configurations  have  limited 
interpretation  of  the  spatially  organized  coherent  energy  to  frequencies  be¬ 
tween  0.  2  and  1.0  Hz. 

Station  power  densities  are  found  to  differ  considerably  across 
the  network  during  each  time  period  and  to  vary  considerably  from  sample  to 
sample.  The  bulk  of  seismic  energy  appears  to  be  concentrated  in  the  0.  2- 
to  0.  5-Hz  frequency  band.  At  1.  0  Hz,  station  power  densities  vary  from 
extremely  low  at  TFO  (-4  to  -8  db)  and  NPNT  (-3  to  -5  db)  to  extremely  high 
at  ADIS  {14  to  26  db)  and  GGGR  (9  to  12  db)  relative  to  1.  0  m|J^/Hz.  These 
variations  in  noise  power  level  with  station  and  time  appear  to  be  directly  re¬ 
lated  to  regional  meteorological  disturbances  and  indicate  a  lack  of  both  time- 
and  space-stationarity  at  the  network  level. 

Station  ambient  noise  fields  during  the  16  October  1964  and  29 
January  1965  time  periods  are  dominated  primarily  by  low- velocity  surface¬ 
mode  energy.  The  exception  is  TFO  where  bodywave-mode  coherent  energy 
is  generally  dominant.  The  coherent  energy  appears  to  be  highly  directional 
and  generally  related  to  regional  storm  activity,  especially  in  the  micro- 
seiemic  band  (2-  to  5-sec  periods).  Wind-generated  surf  action  seems  to  be 
a  primary  source  mechanism.  Storm  activity  is  usual  in  the  northern  latitudes 


1-2 


sclsncs  ssrvlcss  division 


during  the  winter  months,  and  conclusions  from  this  study  may  not  be  appli¬ 
cable  to  other  seasons.  At  frequencies  above  the  microseismic  band,  sig- 
nificant  contribution  from  apparently  stationary  sources  also  is  observed, 
usually  in  the  form  of  localized  narrowband  surface-mode  energy  which  is 
probably  culturally  related. 

While  the  more  intense  storm  sources  generate  energy  simul¬ 
taneously  observed  at  several  stations  in  the  network,  the  amplitude,  fre¬ 
quency,  and  propagation  mode  of  the  observed  energy  vary  with  station.  Such 
variation  is  probably  related  to  several  factors  including  station- to- source 
distances,  source  radiation  patterns,  travel-path  geology  differences,  and 
the  presence  of  strong  regional  sources  that  obscure  weaker  energy  from 
distant  sources.  This  lant  factor  is  particularly  evident  at  the  lower  fre¬ 
quencies  where  array  resolutions  are  poor.  Therefore,  the  low-velocity 
c  ^herent  energy  that  is  a  major  contributor  to  the  noise  field  at  the  station 
level  is  apparently  not  coherent  at  the  network  level  in  the  short-period  fre¬ 
quency  range,  since  t  is  primarily  generated  from  regional  sources.  Wiener 
coherent  processing  does  not  appear  advantageous  for  suppression  of  this 
noise  at  the  network  level. 

Except  for  TFO,  which  previous  studies  show  to  be  bodywave 
noise-limited,  the  array  stations  do  not  appear  to  be  bodywave  noise-limited 
during  the  two  time  periods  of  intense  storm  activity.  However,  bodywave 
noise,  appears  to  form  a  larger  proportion  of  the  total  noise  power  density  at 
NPNT  than  at  CPO,  OONW,  or  GGGR.  Measurements  of  bodywave  noise 
levels  across  the  network  suggest  that  these  levels  are  not  uniform  but  vary 
with  continent.  At  each  station,  the  bodywave  noise  levels  appear  relatively 
inconsistent  with  time  over  the  0.  2-  to  1.  0-Hz  range. 
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Examination  of  the  long-period  data  suggests  valid  seismic 

energy  between  0.  03  and  0.  2  Hz.  A  major  peak  is  seen  around  0.  15  Hz,  with 

a  lesser  peak  indicated  between  0.  05  and  0.  08  Hz.  Spectral  levels  are  notably 

2  , 

consistent  at  0.  04  Hz,  averaging  about  45  db  relative  to  1  mu  / Hz.  No  usable 
coherence  is  observed  between  stations  for  the  long-period  data,  but  this  is 
not  surprising  in  view  of  the  regional  character  of  the  noise  field  already 
reported  for  the  short-period  data;  however,  differences  in  tape  speeds  for 
the  recorded  data  result  in  slightly  different  digitization  rates  at  each  station, 
and  the  coherence  measurements  certainly  are  affected.  Low-level  coher¬ 
ence  (significantly  above  the  level  expected  for  random  noise)  seen  for  cer¬ 
tain  stations  around  0.25  Hz  may  possibly  be  seismic. 

Broadband  spectral  estimates  formed  by  joining  the  absolute 
short-period  and  long-period  power  spectra  appear  generally  reasonable  and 
indicate  a  power  peak  in  the  spectral  range  between  0.  1  and  0.  3  Hz.  Power 
levels  above  0.  2  Hz  appear  to  decay  logarithmically. 

Due  to  the  preponderance  of  surface-mode  energy  at  the  net¬ 
work  stations,  the  extent  to  which  the  bodywave  component  of  the  noise  field 
is  coherent  across  the  network  is  not  established.  With  sufficient  prepro¬ 
cessing  to  reduce  the  low-velocity  component  significantly  at  the  3tation 
level,  useful  coherence  could  conceivably  be  found  for  an  essentially  body- 
wave  noise-limited  network.  In  this  case,  Wiener  coherent  processing 
would  prove  useful  for  suppressing  highly  directional  energy  from  bodywave 
sources. 

Conclusions,  it  should  be  emphasized,  are  drawn  from 
examination  of  a  very  limited  quantity  of  data;  results  cannot  be  generalized 
for  year- round  network  operation  without  further  analysis  using  data  samples 
representative  of  all  seasons.  Also,  the  insufficient  number  of  network  sta¬ 
tions  available  precluded  comprehensive  coverage.  Continued  evaluation  with 
an  expanded  network  of  quality  array  stations  and  an  ensemble  of  noise  sam¬ 
ples  representing  year-round  network  operation  is  recommended. 
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SECTION  II 


NETWORK  NOISE  POWER  SPECTRAL  ANALYSIS 
A.  SHORT -PERIOD  NOISE  POWER  SPECTRA 

Power-density  spectra  are  computed  using  available  noise 
samples  at  each  station.  Single  vertical- component  short-period  seismom¬ 
eters  are  used  to  represent  each  station.  These  spectra  are  then  used  to 
interpret  data  validity  as  a  function  of  frequency,  the  time-  and  space- 
stationarity  of  network  noise,  location  of  spectral  peaks,  and  general  dis¬ 
tribution  of  power  as  a  function  of  frequency. 

Data  gate  lengths  vary  from  370  sec  to  1400  sec,  with  data 
quality  a  determinir"  factor.  Trace  areas  containing  spikes  or  other  ob¬ 
viously  nonseismic  disturbances  have  been  avoided.  Autocorrelations  are 
computed  for  ±10  sec.  The  correlations  are  Parzen-smoothed,  yielding 
a  frequency  resolution  of  approximately  0.  1  Hz  in  the  resulting  power  spectra. 

Figures  II- 1  through  II- 9  present  station  power  spectra  (cor¬ 
rected  only  at  1.0  Hz)  for  the  various  network  stations.  Available  spectra 
at  each  station  are  overlaid,  facilitating  comparison  of  spectral  variation 

with  time.  The  power  variation  in  the  microseismic  band  with  season  shows 

3 

general  agreement  with  previous  work.  Station  noise  levels  appear  to  vary 
primarily  as  a  function  of  regional  storm  activity,  and  a  general  rise  in 
noise  power  is  evident  across  the  network  during  the  winter  months.  Since 
the  network  consists  of  northern-hemisphere  stations  only,  this  rise  is  in¬ 
dicative  of  increased  storm  activity  in  the  northern  hemisphere  during  the 
winter  months. 
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Noise  energy  levels  at  ADIS  appear  to  vary  considerably  with 
time  (Figure  II- 1).  Most  notable  is  the  considerably  reduced  level  on  12  Aug¬ 
ust  1964.  The  elevated  level  on  6  November  1964  reflects  the  presence  of 
an  intense  low-pressure  center  (LPC)  off  the  Aleutian  chain.  Background 
levels  are  generally  higher  than  for  stations  located  on  the  continent,  sug¬ 
gesting  that  surf  action  may  be  a  primary  source. 

Figure  II- 2  shows  considerable  energy  at  CPO  for  0.  2  Hz 

during  16  October  1964,  indicating  the  presence  of  tropical  storm  Isabel 

along  the  Atlantic  coast.  The  spectra  at  CPO  appear  generally  well- equalized. 

Noise  power  appears  slightly  above  0  db  at  1.  0  Hz,  agreeing  well  with  pre- 
3 

vious  work.  There  are  evidences  also  of  spectral  peaks  at  1.4  and  2.  0  Hz; 
both  correspond  to  previously  reported  directional  surface-mode  energy. 

High- resolution  spectral  analysis  suggests  that  these  peaks  correspond  to 
a  common  local  source  north-northwest  of  CPO. 

An  interesting  characteristic  of  the  spectra  at  GGGR  is  evi¬ 
dent  in  Figure  II- 3 .  While  microseismic  power  levels  are  generally  low 
relative  to  the  other  stations,  the  1.4-  to  2.  5-Hz  frequency  band  reflects 
unusually  high  power  levels  with  numerous  peaks;  whether  these  peaks  are 
seismic  is  not  established. 

Figure  II-4  shows  the  power  spectra  at  OONW.  The  spectrum 
for  6  November  1964  indicates  a  much  higher  overall  noise  level  up  to  2.  0 
Hz.  An  intense  low-pressure  center  was  located  over  North  Europe  during 
this  time  period.  The  12  August  1964  spectrum  is  again  relatively  lower 
across  the  microseismic  band,  indicating  a  lack  of  regional  storm  activity. 

All  spectra  are  from  compressed  FM  data,  except  the  16  October  1964  sam¬ 
ple,  which  exhibits  generally  lower  levels  above  2.  0  Hz. 
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Figure  II- 1.  Short-Period  Vertical-Component  Power  Spectra  at  ADIS 


Figure  II-2.  Short-Period  Vertical- Component  Power  Spectra  at  CPO _ 
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Figures  II- 5  and II - 6  for  RKON  and  DHNY,  respectively,  indi¬ 
cate  generally  lower  microseismic  energy  levels  on  12  August  1964.  Figures 
II- 7 ,  II- 8 ,  and  II- 9  contain  spectra  computed  for  NPNT,  LZBV,  and  TFO, 
respectively.  The  16  October  1964  spectrum  at  TFO  appears  odd  and  is  not 
believed  to  be  seismically  valid. 

In  general,  the  16  October  1964  and  29  January  1965  samples 
are  from  noncompressed  data  (except  OONW  on  29  January  1965)  and  reflect 
lower  levels  above  2.  0  Hz.  All  others  are  from  compressed  recordings  and 
indicate  numerous  spectral  "lines"  as  well  as  a  generally  higher  power  level 
at  frequencies  above  2.  0  Hz. 

Figure  II- 10  shows  the  power-density  level  at  1.  0  Hz  eeen  at 
each  station  for  several  noise  samples.  TFO  reflects  the  lowest  levels,  fol¬ 
lowed  closely  by  NPNT  and  DHNY.  CPO,  LZBV,  RKON,  and  OONW  appear 
intermediate;  GGGR  and  ADIS  indicate  relatively  high  ambient  noise  levels. 
The  6  November  1964  sample  is  consistently  higher  i:  noise  power  across 
the  network;  other  samples  indicate  more  variance. 

Calibration  analysis  performed  on  available  calibration  data 
is  used,  together  with  station-log  data,  to  determine  correction  factors  sub¬ 
sequently  applied  to  the  power  spectra,  yielding  power  density  in  decibels 
2 

relative  to  1  m|i  /Hz  at  1.  0  Hz. 

This  procedure  encountered  several  problems,  however: 
station  logs  were  not  available  for  all  samples,  and  the  available  ones  were 
sometimes  ambiguous;  in  a  few  instances,  calibrations  were  not  available  or 
usable  for  analysis.  Estimates  from  other  days  were  used  in  several  cases. 

Some  corrections,  although  consistent  with  available  infor¬ 
mation,  appear  to  provide  questionable  normalization  of  spectra.  In  general, 
the  overall  quality  of  the  normalization  through  the  use  of  calibrations  is  less 
than  desired. 
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Several  conclusions  can  be  drawn  from  examination  of  the 
spectra,  however.  All  stations  except  GGGR  appear  dominated  by  the  0.2- 
to  0.  5 -Hz  microseismic  energy.  Analysis  using  f-k  spectra  (discussed  in 
Sections  III  and  IV)  indicates  that  the  bulk  of  this  energy  is  of  low  velocity 
(surface  mode);  the  exception  is  TFO.  GGGR  appears  somewhat  unique  in 
that  spectral  noise  power  is  significantly  higher  in  the  signal  band,  but  no 
explanation  for  this  elevated  power  level  is  available. 

•*» 

Above  2.  0  Hz,  the  spectra  rapidly  break  into  narrowbai  d 
peaks.  Most  of  these  peaks  are  not  believed  to  be  seismic  but  rather  to  be 
generated  during  FM  recording  or  playback,  during  compression  of  data, 
or  during  digitization.  Extremely  high-velocity  energy  is  indicated  by  f-k 
analysis  for  these  peaks.  Ideally,  instrument-generated  noise  might  be 
expected  to  exhibit  infinite-velocity  peaks.  However,  the  process  of  digiti¬ 
zation  introduces  a  small  but  finite  time  delay  in  sampling  between  channels, 
which  appears  in  f-k  analysis  as  a  relatively  consistent  phase  shift  across 
the  array.  This  timing  error  is  approximately  1/600  sec  between  adjacent 
channels  for  standard  FM-to-digital  tapes  and  1/120  sec  for  compressed 
FM-to-digital  tapes. 

The  strong,  very  narrowband  peak  at  3.6  Hz  seen  in  the 
OONW  6  November  1964  spectrum  has  been  analyzed  through  direct  trans¬ 
forms  computed  for  all  array  channels.  Figure  II- 11  shows  the  phase  shifts 
observed  for  each  channel  relative  to  Z2.  The  straight  line  represents 
theoretical  phase  shifts  expected  from  compressed  data,  assuming  a  1/ 120- 
sec  delay  in  sampling  adjacent  channels  containing  in-phase  information  at 
3.6  Hz.  Except  for  Zl,  which  is  lowest  in  amplitude  at  3.6  Hz,  there  is 
excellent  agreement  between  measured  and  theoretical  phase  shifts.  There¬ 
fore,  this  timing  error  inherent  in  sampling  is  interpreted  as  the  mechanism 
for  producing  the  deviations  from  infinite  velocity  seen  for  instrumentally 
produced  spectral  peaks. 
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SEIS  NO. 

Figure  II- 11,  Comparison  between  Measured  and  Theoretical 
Phase  Shifts  Relative  to  Channel  Z2,  OONW 

The  spectra  generally  are  not  believed  seismically  valid  above 
2.0  Hz.  Similarly,  estimates  below  0.2  Hz  are  questionable.  Therefore, 
the  frequency  range  between  0.2  and  2.0  Hz  appears  to  be  generally  useful 
for  network  spectral  analysis. 

Based  on  the  limitad  amount  of  data  analyzed,  network  power 
spectra  appear  neither  time-  nor  space-stationary.  The  apparent  dependence 
of  station  power  level  on  regional  effects  suggests  strongly  an  overall  lack  of 
coherence  at  the  network  level.  Further  study  is  indicated,  using  analysis  of 
a  wide  range  in  data  samples  representing  various  seasonal  and  diurnal  con¬ 
ditions.  If  the  frequency  range  of  valid  seismic  data  is  to  be  extended,  higher- 
quality  network  data  probably  will  be  required. 


B.  LONG- PERIOD  SPECTRAL  ANALYSIS 


Figure  II- 12  shows  long-  period  power-density  spectra  com¬ 
puted  from  vertical- component  data  for  16  October  1964  and  29  January  1965. 
Pa  rzen- smoothed,  the  correlations  used  in  determining  the  spectra  are  com¬ 
puted  for  100  lags  (50  sec),  yielding  a  frequency  resolution  of  approximately 
0.  02  Hz.  Simultaneous  data  gates  approximately  30  min  in  length  are  used 
from  each  station. 

The  power  spectra  are  shown  corrected  to  units  of  absolute 
ground  motion  at  0.  04  Hz  for  those  samples  where  calibration  information 
is  available.  Also  shown  in  subsection  C  are  several  spectra  with  instru¬ 
ment  responses  removed  for  the  broadband  spectra  estimates.  The  spectra 
at  TFO  and  HWIS  are  shown  relative  to  unity,  since  calibration  information 
is  not  available  for  these  samples. 

Three  major  peaks  characterize  the  power  spectra  for  both 
noise  samples.  The  sharp  peak  at  0.6  Hz  is  nonseismic  and  probably  was 
introduced  during  FM  data  compression  operations  (or  FM-to-digital  con¬ 
version).  However,  the  stronger  peaks  at  0.  04  to  0.  06  Hz  and  at  0.  12  to 

0.  14  Hz  do  appear  seismic  and  in  close  agreement  to  those  reported  at 
4 

LAS/' .  A  third  peak,  which  occurred  near  0.  19  to  0.21  Hz  at  CPO,  ADIS, 

and  Ha'IS  on  16  October  1964,  corresponds  to  the  broad  peak  from  0.  12  to 

0.22  Hz  observed  at  LASA.  The  spectra  appear  to  contain  generally  valid 

seitmic  data  between  0.  03  and  0.  2  Hz  and  seem  to  be  very  stationary  at 

0°  04  Hz,  with  all  stations  reflecting  levels  very  close  to  45  db  relative 
2 

to  1  m^i  /Hz.  Larger  variations  are  seen  for  the  microseismic  energy 
between  0.  1  and  0.  2  Hz.  The  rise  in  power  below  0.  03  Hz  is  believed  re¬ 
lated  to  instrument  or  meteorological  energy  rather  than  to  seismic  energy. 
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In  general,  the  quality  of  the  long-period  data  is  low.  Major 
portions  of  the  recorded  data  are  unusable  because  of  spikes.  Data  digiti¬ 
zation  levels  are  very  low,  causing  rather  coarsely  quantized  samples  and 
poor  seismic-to-instrument-noise  ratios.  A  third  problem  affecting  coher¬ 
ence  measurements  involves  apparent  tape- speed  differences  in  the  station- 
recorded  data  during  digitization  of  the  FM  data  (or  during  original  recording 
of  the  FM  tapes).  These  differences  result  in  slightly  different  sampling 
rates  for  each  station.  For  example,  sampling  rates  on  29  January  1965 
vary  from  7620  samples /hr  at  ADIS  to  8100  samples /hr  at  CPO,  based  on 
1  hr  of  station- recorded  clock  time;  this  contrasts  with  the  "assumed"  7200 
samples/hr  (2  samples/sec).  Total  sampling-rate  error  is  higher  for  the 
16  October  1964  data,  but  station-to- station  variations  are  smaller  (e.  g. , 
8600  samples /hr  at  RKON  in  contrast  with  8500  samples /hr  at  CPO). 

Below  0.06  Hz,  observations  are  necessarily  tentative,  since 
data  gates  were  only  30  min  long.  However,  the  apparently  stable  character 
of  the  spectra  at  0.  04  Hz  implies  a  relative  lack  of  sensitivity  to  regional 
seismic-noise  activity.  Additional  evaluation  of  this  behavior  is  indicated. 

Figure  11-13  shows  interstation  coherences  computed  for 
the  stations  included  in  the  power  spectral  analysis.  Except  for  highly 
coherent  lines  at  0.6  Hz  and  broad  peaks  below  0.  03  Hz  (all  nonseismic), 
usable  coherence  is  nil  across  the  network.  The  sampling  error  described 
earlier  certainly  affects  these  measurements;  however,  in  view  of  the  dom¬ 
inance  of  regional  surface-mode  energy  at  each  station  in  the  short-period 
range,  little  coherence  would  be  expected  in  the  microseismic  band. 

Assuming  that  the  long- period  network  data  correspond  to 
a  weakly  st?tionary  process,  an  estimate  of  the  level  of  coherence  expected 
through  CiMtnce  correlation  can  be  made.  From  tables  of  the  distribution 
of  coherence  coefficients  presented  by  Amos  and  Koopmans,  5  this  level  is 
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0.  15,  with  a  standard  deviation  of  0.  075.  The  coherence  curves  presented 
in  Figure  11-13  actually  represent  coherence  squared  and  are  computed  using 
the  relation 
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where  *  denotes  the  complex  conjugate.  Since  tables  for  the  distribution  of 
estimates  for  coherence  squared  are  not  as  readily  available,  values  of  in¬ 
terest  are  converted  to  Koopmans'  definition  by  taking  the  square  root. 

Several  coherent  peaks  are  significantly  higher  than  the 
random  correlation  level;  whether  these  are  instrument-generated  or  seis¬ 
mic  is  difficult  to  establish.  Especially  interesting,  however,  is  one  peak 
representing  0.4  coherence  (0.  15  on  plot  of  coherence  squared)  at  0.26  Hz 
between  RKON  and  TFO.  Both  spectra  are  substantially  above  the  background 
digitization  levels  (Figure  11-12)  at  this  frequency  and  both  show  slight  evi¬ 
dence  of  a  spectral  peak.  Also,  the  two  stations  are  essentially  in  line  with 
coherent  energy  propagating  across  the  Canadian  Shield  from  the  intense 
storm  located  off  Newfoundland.  Such  coherence,  if  seismic,  is  too  low  for 
practical  use  in  array  processing  but  certainly  is  academically  interesting. 
Also,  this  coherence,  if  truly  seismic,  probably  would  be  increased  by  the 
analysis  of  larger  time  samples. 

In  summary,  the  long-period  spectra  appear  to  be  seismically 
valid,  to  reflect  reasonably  stable  energy  levels  in  the  0.  03-  to  0.  05-Hz  band 
where  the  instrument  responses  peak,  and  to  indicate  that  microseismic  energy 
is  a  significant  contributor  between  0.  1  and  0.  3  Hz.  Energy  levels  in  this  lat¬ 
ter  band  appear  to  be  related  to  regional  storm  activity.  No  usable  coherence 
between  network  stations  is  detected;  however,  due  to  the  limited  quality  of 
the  data,  conclusions  regarding  network  coherence  in  the  long-period  range 
cannot  be  drawn. 
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Figure  11-14.  Broadband  Spectral  Estimates  Using  Absolute 
Long-  and  Short-Period  Power-Density  Spectra 
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C.  BROADBAND  SPECTRAL  ESTIMATES 
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Figure  11-14  presents  absolute  power-density  spectra  for  CPO, 
NPNT,  ADIS,  and  RKON.  The  frequency  range  from  0.  01  to  2.0  Hz  is  shown, 
with  long-period  spectral  estimates  used  to  form  the  curves  from  0.  01  to  0.2 
Hz  and  short-period  data  used  from  0.  2  to  2.  0  Hz. 

The  spectra  join  reasonably  well  in  the  0.  2-Hz  vicinity  and 
appear  to  provide  reasonable  estimates  of  the  broadband  seismic  noise  field. 
The  seismic  peak  seen  on  all  spectra  between  0.  1  and  0.  3  Hz  is  apparently 
valid,  although  the  exact  level  is  suspect  since  both  long-period  and  short- 
period  instrument  responses  are  low  (approximately  -40  db)  at  this  point 
and  seismic-noise-to-instrument-noise  levels  are  poor.  Below  0.03  Hz, 
the  long-period  curves  are  also  suspect  due  to  contamination  with  system 
noise.  Noise  power  levels  are  up  considerably  at  0.  04  Hz  (45  db)  relative 
to  1.0  Hz  (-10  to  +10  db)  and  are  much  stabler  across  the  network. 
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SECTION  III 


ANALYSIS  OF  NETWORK  NOISE  FIELD  USING  HIGH- 
RESOLUTION  FREQUENCY-WAVENUMBER  SPECTRA 

A.  DESCRIPTION  OF  EXPERIMENT 

A  knowledge  of  the  noise  structure  at  the  network  level  is 
necessary  for  optimum  design  of  signal  detection  and  extraction  procedures 
for  the  network.  Therefore,  an  objective  of  the  network  noise  study  is  dis¬ 
section  of  each  station’s  short-period  noise  field  into  frequency  and  wave- 
number  components  and  determination  of  component  sources  where  possible. 
The  network  noise  structure  may  then  be  interpreted  in  terms  of  energy  from 
common  sources  when  seen  at  each  station.  The  network's  capability  for  lo¬ 
cating  and  tracking  energy  sources  such  as  tropical  storms  can  also  be  evalu¬ 
ated. 


Such  an  analysis  requires  reasonably  long  simultaneous  noise 
samples  from  each  array  station  in  the  network.  Practical  limitations,  in¬ 
cluding  data  quality,  resulted  in  a  usable  data  gate  length  of  6.  67  min  at  each 
station.  Two  of  the  five  available  noise  samples  were  taken  during  periods 
when  low-pressure  areas  were  numerous  across  the  northern  hemisphere; 
these  samples  (16  October  1964  and  29  January  1965)  were  chosen  for  anal¬ 
ysis.  Data  for  these  samples  were  available  from  array  stations  CPO,  TFO, 
OONW,  GGGR,  and  NPNT .  Descriptive  material  pertaining  to  these  stations 
and  the  data  library  is  contained  in  the  appendix. 

A  technique  useful  for  dissection  must  provide  good  resolution 
in  both  frequency  and  wavenumber.  Therefore,  the  high- resolution  f-k  spec¬ 
tral  technique  was  employed  in  conjunction  with  direct-transform  procedures. 
The  use  of  these  techniques  for  noise  analysis  was  evaluated  in  a  previous 
special  report.  ^ 

Smoothing  across  frequency  was  not  performed  on  the  resulting 
crosspower  matrices,  so  frequency  resolution  was  roughly  0.0025  Hz.  Com¬ 
parisons  described  in  Special  Report  No.  8  demonstrate  that  phase  relation¬ 
ships  between  sensors  are  distorted  by  smoothing  over  long  time  gates,  re- 
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suiting  in  increased  variability  in  spatial  location  of  peaks  with  choice  of 
reference  sensor.  Further  comparisons  of  the  unsmoothed  spectra  with  the 
rather  well-documented  CPO  noise  field  (in  terms  of  conventional  f-k  anal- 

•j 

ysis)  indicate  excellent  agreement. 

Since  the  high- resolution  technique  is  influenced  by  reference- 
sensor  choice  (due  to  the  finite  array  geometry),  wavenumber  responses  were 
computed  using  each  array  sensor  as  reference.  These  MCF  responses  are 
actually  the  reciprocal  of  the  wavenumber  spectral  estimate.  Prior  to  the 
summing  across  sensors,  the  reciprocal  of  each  response  was  taken,  yield¬ 
ing  additional  discrimination  between  coherent  peaks  and  the  overall  incoher¬ 
ent  level  of  the  plane.  Crosspower  matrices  were  normalized  to  yield  unit 
autopower  at  each  frequency.  A  signal-to-noise  ratio  of  0.  01  used  in  the  fil¬ 
ter  design  was  experimentally  determined  adequate  for  resolution  in  the  wave- 
number  plane. 

An  incentive  for  critical  analysis  of  high- resolution  spectra 
computed  from  unsmoothed  direct  transforms  taken  over  rather  long  time 
gates  is  the  emergence  of  a  mathematical  shortcut,  facilitating  computation 
of  the  high- resolution  spectral  sum.  The  procedure  is  valid  only  for  un¬ 
smoothed  crosspower  matrices  and  eliminates  the  need  to  compute  and  sum 
individual  filter  responses  for  each  sensor.  Using  this  shortcut,  the  summed 
output  was  obtained  in  roughly  the  time  originally  required  for  each  individual 
spectrum.  Special  Report  No.  8  contains  a  description  of  the  mathematical 
procedure  used.^ 

As  a  result  of  summing  reciprocals  of  the  individual  wave- 
number  responses,  some  smoothing  in  the  wavenumber  plane  appears  de¬ 
sirable.  This  is  caused  primarily  by  departures  from  ideal  equalization 
of  coherent  energy,  especially  at  the  lower  frequencies,  and  is  seen  in  the 
individual  spectra  as  slight  wandering  of  the  coherent  peaks  in  the  wavenum¬ 
ber  plane.  Procedures  to  smooth  in  the  wavenumber  plane  require  further 
examination  and  were  not  implemented  for  this  study  Instead,  a  minimum 
of  hand- smoothing  was  performed  during  contouring,  primarily  to  enhance 
th*  spectral  readability. 
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High- resolution  f-k  spectra  were  computed  for  each  station 
in  0.2-Hz  increments  from  0.  2  Hz  to  1 . 0  Hz  (TFO),  1.4  Hz  (NPNT,  OONW, 
GGGR),  and  2.  0  Hz  (CPO).  Above  2.  0  Hz,  the  seismic  validity  of  the  net¬ 
work  data  is  questionable,  especially  with  the  compressed  FM  data.  Also, 
minimum  seismometer  spacings  at  the  array  sites  have  generally  led  to 
space-aliasing  of  valid  low-velocity  peaks  at  the  higher  frequencies,  making 
interpretation  impractical. 

B.  RESULTS  AND  CONCLUSIONS 

The  spatially  organized  noise  observed  at  the  station  level 
appears  to  be  generated  primarily  by  meteorological  and  cultural  sources. 

The  former  source  type  includes  regional  weather  effects  (such  as  nearby 
low-pressure  centers,  frontal  passages,  and  wind-generated  surf  action) 
and  produces  the  bulk  of  the  coherent  energy.  The  cultural  sources  appear 
stationary  and  highly  localized,  generating  low-velocity  (surface-mode) 
energy  that  appears  random  at  the  network  level.  Above  the  microseismic 
range,  these  local  sources  are  substantial  contributors  to  the  station  noise. 

Weather-generated  energy  is  a  continuous  contributor  through¬ 
out  the  frequency  band  (to  2.  0  Hz)  and  dominates  in  the  microseismic  range. 
Table  III- 1  summarizes  coherent  energy  apparently  related  to  regional  storms; 
subsection  C  presents  and  discusses  the  station  high- resolution  f-k  spectra 
containing  these  coherent  peaks. 

Major  coherent  peaks  at  TFO  are  primarily  bodywave  energy; 
CPO  and  the  LRSM  stations  are  dominated  by  low-velocity  energy.  Contri¬ 
buting  to  the  sensitivity  at  TFO  are  several  factors  including  array  aperture 
(10  km  as  compared  with  3  km  at  CPO  and  the  LRSM  stations),  geology  (scat¬ 
tering  of  surface  waves  at  TFO  has  been  noted),  relative  ambient  noise  levels, 
and  a  lack  of  localized  storm  activity  in  the  TFO  vicinity. 


Ill- 3 


solsnes  ssrvlcss  division 


Table  III- 1 

STORM-GENERATED  COHERENT  ENERGY 
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Nearby  storm  centers  (Buch  as  LPC  1  east  of  CPO  and  LPC  13 
in  the  North  Sea  west  of  OONW)  radiate  energy  over  a  broad  frequency  band. 
However,  the  coherent  peaks  generally  appear  rather  narrowband.  This  is 
due  partly  to  the  high-frequency  resolution  available  —  negligible  energy  is 
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smeared  across  frequencies.  In  addition,  the  attenuation  of  energy  with  both 
distance  and  geology  is  frequency-dependent.  The  weather- related  seismic- 
energy  source  functions  are  not  necessarily  broadband. 

While  the  more  intense  low-pressure  centers  are  simultaneously 
detected  at  several  network  stations,  this  detection  is  rarely  achieved  at  a  common 
frequency  and  only  once  for  bodywave  energy  (LPC  20  seen  at  both  OONW  and  GGGR 
for  1 . 2  Hz).  Simultaneous  detection  of  surface- wave  energy  occurs  for  LPC  20  (CPO 
*ud  .'NT  at  2.  0  Hz),  LPC  22  (NPNT  and  GGGR  at  0.  4  and  1 .  0  Hz),  LPC  14  (NPNT  and 
TFOat  0.4  and  1.0  Hz),  and  LPC  13  (OONW  and  GGGR  at  0.  2  and  1 . 0  Hz).  This  lack  of 
common  detection  is  partly  explained  by  the  attenuation  of  energy  with  distance  to  each 

station  as  a  function  of  frequency.  Also,  differences  infrequency  filtering  with  azi- 
muth  at  each  station  contribute  to  this  effect.  Similarly,  source- radiation  patterns 
probably  are  not  uniform  with  azimuth. 


The  CPO  noise  spectra  are  compared  with  data  from  a  rather 
comprehensive  noise  analysis  using  conventional  f-k  spectra.7  This  is  done 
partly  to  confirm  the  validity  of  the  data  and  partly  to  examine  the  validity  of 
the  high- resolution  f-£  technique  for  such  noise  analysis.  Excellent  agree¬ 
ment  is  obtained  in  terms  of  fe^ral  spatial  organization  and  specifically  for 
established  stationary  sourc,  .  >.g.  ,  the  narrowband  coherent  peak  seen  at 
1.4  Hz  from  the  north  at  about  4.  5  km/sec).  Similarly,  general  spatial  dis¬ 
tribution  at  TFO  agrees  well  with  previous  analyses.  8 


A  few  high-velocity  peaks  which  do  not  appear  to  be  either 
storm-generated  or  stationary  are  detected.  Sources  for  these  peaks  are 
not  established,  although  there  is  the  possibility  that  they  represent  valid 
bodywave  energy  from  seismically  active  areas. 


The  foregoing  analysis,  it  must  be  emphasized,  is  based  on 
a  rather  limited  quantity  of  network  data.  Only  two  noise  samples,  both 
highly  storm-active,  are  used.  More  realistically,  several  samples  cover¬ 
ing  a  range  of  weather  and  seasonal  conditions  should  be  analyzed.  The  net¬ 
work  of  five  array  stations  used  actually  is  not  sufficient  for  comprehensive 
coverage  and  should  be  expanded. 


Ill- 5 


sclsncs  ssrvleea  division 


Some  preliminary  conclusions  regarding  the  character  of  the 
network  noise  field  can  be  made,  however,  from  the  foregoing  analysis. 

The  lack  of  common  coherent  peaks  at  the  network  stations 
for  a  given  frequency  indicates  that  little  or  no  coherence  should  be  found 
at  the  network  level.  Thus,  there  appears  to  be  little  advantage  in  using 
multichannel  processing  techniques  for  noise  suppression  at  the  network 
level.  The  bulk  of  the  station  noise  appears  as  low-velocity  energy  (except 
at  TFO),  and  no  conclusions  are  drawn  regarding  the  extent  of  coherence  of 
the  bodywave  portion  of  the  station  noise  field.  Generally,  the  station  arrays 
are  not  adequate  to  sufficiently  suppress  the  low-velocity  component,  allowing 
evaluation  of  bodywave  coherence  across  the  network.  Should  preprocessing 
at  the  station  level  provide  suppression  of  the  low-velocity  component,  multi¬ 
channel  processing  might  conceivably  be  required  for  directional  bodywave 
noise  suppression  at  the  network  level. 

C.  PRESENTATION  OF  STATION  NOISE  SPECTRA 

High- resolution  spretra  from  the  five  array  stations  are  ana¬ 
lyzed  at  the  same  time  for  each  frequency.  Data  from  the  16  October  1964 
noise  sample  are  presented  first,  followed  by  the  29  January  1965  data. 

Since  a  rather  large  number  of  spectra  are  presented,  the  following  format 
facilitates  discussion  of  both  individual  and  common  features:  spectra  from 
the  five  stations  are  shown  together  at  each  frequency,  and tbediscussicn  per¬ 
tinent  to  these  spectra  is  contained  on  the  opposite  page.  Shown  with  each 
set  of  spectra  is  a  map  displaying  the  relative  locations  of  stations  and  storm 
centers  for  that  day. 

The  wavenumber  plaue  out  to  an  edge  velocity  of  2.  0  km/sec 
is  shown.  At  the  lower  frequencies,  this  region  is  entirely  contained  within 
the  unit  cell.  However,  as  frequency  increases,  the  available  wavenumber 
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region  is  reduced  by  contraction  of  the  unit  cell.  Dimensions  of  each  unit 
cell  are  indicated  when  they  fall  within  *he  2. 0-km/sec  region  displayed. 

The  plots  are  not  contoured  outside  the  unit  cell  except  for  a  few  instances 
when  a  low-velocity  coherent  peak  is  believed  aliased  into  the  cell  from 
outside.  Differences  in  minimum  eeismometer  spacings  along  the  two  arms 
at  TFO  result  in  a  rectangular  unit  cell. 

Contour  levels  vary  from  plot  to  plot  due  to  variations  in 
dynamic  range  and  complexity.  Within  each  plot,  however,  the  contour  in¬ 
crement  is  held  constant.  Exception'!  occur  when  relative  peak  amplitudes 
are  indicated.  All  plots  are  oriented  conventionally,  with  north  at  the  top 
of  the  page  (except  CPO,  where  true  north  lies  14°  west  of  conventional 
north)  as  shown  by  the  arrow  at  the  top  of  each  plot.  Coherent  spectral 
peaks  occur  in  the  quadrant  corresponding  to  the  direction  of  the  source 
from  the  station.  Shading  is  used  to  highlight  peaks  being  discussed.  Also, 
the  peaks  are  assigned  numbers  corresponding  to  storm-center  designations 
if  believed  to  be  related  to  that  source.  Apparently  stationary  peaks  are 
designated  with  letters. 

The  discussion  refers  to  storm  centers  as  low-pressure  cen¬ 
ters  (LPC).  Thus,  tropical  storm  Isabel,  designated  storm  center  1,  is 
referred  to  as  LPC  1.  Figures  III- 1  and  HI- 2  show  areas  of  storm  activity 
for  the  northern  hemisphere  on  16  October  1964  at  1800  GMT  and  29  January 
1  ?65  at  0600  GMT,  respectively.  Investigation  of  southern-hemisphere  wea¬ 
ther  data  indicates  no  significant  storm  activity.  Isobaric  lines  of  high-  and 
low-pressure  areas  and  frontal  lines  are  shown  on  each  weather  map.  Fig¬ 
ures  III- 3  and  IU-4  show  sea  wind-wave  height  contours  and  indicate  the  lo¬ 
cation  of  designated  LPC's.  Swell-height  charts  reflect  similar  ^ea-state 
information  and  are  not  shown.  Figures  III- 5  t  trough  III-  19  contain  the  2- 
dimensional  wavenumber  spectra. 
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Figure  III-3.  Sea  Wind-Wave  Height  Chart  for  16  October  1964  at  1800  GMT 
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Figure  HI-4.  Sea  Wind- Wave  Height  Chart  for  29  January  1965  at  0600  GMT 
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Figure  III-5,  16  October  1964,  F  =  0.  2  Hz 

Surface- mode  energy  generated  by  regional  storm 
activity  dominates  the  station  noise  fields  at  0.  2  Hz  during  this  time 
period.  An  exception  is  TFO,  which  indicates  bodywave  energy  (15  km/ 
sec)  arriving  from  the  northeast  (lobe  14),  apparently  generated  by  the 
intense  LPC  14  located  off  the  Newfoundland  coast. 

The  surface-mode  peak  at  CPO  (lobe  1 ) is  generated  by 
LPC  1  near  Cape  Hatteras.  This  energy  is  observed  over  a  relatively 
wide  frequency  band,  probablv  due  to  the  close  proximity  of  the  distur¬ 
bance. 

Surface-mode  energy  indicated  at  NPNT  (lobe  3)  is  ap¬ 
parently  generated  by  LPC  3  off  the  southern  Alaskan  coast.  GGGR  and 
OONW  indicate  surface-mode  energy  (lobe  13)  from  LPC  13  located  in 
the  North  Sea  vicinity. 

Surface-mode  velocities  vary  from  2.0  to  3.0  km/sec. 
Spectral-peak  definition  is  generally  poor,  primarily  because  of  insuf¬ 
ficient  array  resolution  at  this  frequency. 
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Figure  III-6,  16  October  1964,  F  =  0.4  Hz 


Regional  storms  again  dominate  most  spectra  at  0.4  Hz, 
although  some  seemingly  stationary  peaks  emerge. 

Bodywave  energy  is  seen  at  OONW  (lobe  13),  with  the 
8  km/sec  peak  apparently  generated  by  LPC  13.  However,  evidence 
of  this  source  at  GGGR  is  submerged  by  the  strong  surface-mode  energy 
from  the  south  (lobe  A),  which  corresponds  to  the  major  peak  in  the  sta¬ 
tion  power  spectra  (Section  II).  This  peak  occurs  also  on  29  January 
1965  and  is  probably  related  to  local  cultural  activity.  The  town  of 
Nuremburg,  Germany,  is  28  km  a  vay  in  this  general  direction. 

Surface- mode  energy  from  LPC  1  still  dominates  the 
spectrum  at  CPO  (lobe  1),  and  bodywave  energy  from  LPC  14  continues 
at  TFO  (lobe  14).  Both  sources  are  dominant  contributors  at  0.2  Hz. 
Surface-mode  energy  at  NPNT  (lobe  14)  also  apparently  is  related  to 
LPC  14. 

Coherent  peaks  begin  to  appear  relatively  narrowband, 
suggesting  the  possibility  of  highly  frequency-selective  earth  filtering. 
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Figure  III-7,  16  October  1964,  F  =  0.  6  Hz 


At  0.6  Hz,  storm- related  energy  appears  generally  less 
dominant.  Coherent  energy  seems  relatively  low  at  all  stations,  sug¬ 
gesting  a  transition  region  in  the  0.  6-  to  0.  8-Hz  interval.  For  example, 
coherent  energy  from  LPC  1  has  faded  at  CPO;  a  relatively  broad,  seem¬ 
ingly  stationary  peak  from  the  southeast  (lobe  B)  appears.  The  relatively 
high  surface-mode  propagation  velocity  suggests  a  local  cultural  source. 
This  peak  is  slightly  better  defined  in  the  29  January  1965  sample. 

The  bodywave  peak  seen  at  OONW  at  0.  4  Hz  is  replaced 
by  a  surface-mode  peak  (lobe  13),  probably  generated  by  LPC  13.  Sim¬ 
ilarly,  the  narrowband  stationary  surface-mode  peak  seen  at  GGGR  at 
0.  4  Hz  is  replaced  by  a  peak  from  the  southwest  (lobe  C),  which  also 
occurs  on  29  January  1965  and  is  probably  locally  or  coastally  generated. 
A  weaker  peak  from  the  northeast  (lobe  D)  also  may  be  stationary. 

Bodywave  energy  is  seen  at  NPNT  (lobe  5),  corresponding 
in  wavenumber  to  the  intense  LPC  5  off  Japan.  The  bodywave  peak  at 
TFO  (20  km/sec)  from  the  southwest  (lobe  E)  does  not  appear  to  be  storm¬ 
generated,  and  its  source  is  not  established.  Surface-mode  energy  from 
the  northeast  (lobe  14)  is  apparently  related  to  LPC  14.  A  surface-mode 
oeak  from  the  southeast  (lobe  F)  is  generated  either  locally  or  along  the 
Gulf  coast. 
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At  0.8  Hz,  surface-mode  energy  dominates  at  all  stations. 
Coherent  energy  levels  are  still  relatively  low  at  this  frequency.  GGGR 
indicates  a  rurface-mode  peak  (lobe  13)  from  LPC  13  in  the  North  Sea. 
This  energy  is  indicated  at  0.  2  Hz  and  is  then  replaced  by  stronger  sta¬ 
tionary  peaks  at  0.  4  and  0.  6  Hz. 

Also  seen  is  a  weaker  peak  from  the  south  (lobe  0),  ap¬ 
parently  neither  stationary  nor  storm- related.  OONW  indicates  an 
apparently  stationary  peak  from  the  southeast  (lobe  M).  Lesser  peaks 
from  the  northwest  (lobe  12)  and  south  (lobe  13)  likely  correspond  to 
surf  action  from  LPC's  12  and  13,  respectively. 

The  CPO  spectrum  closely  resembles  that  seen  for 
29  January  1965,  with  several  apparently  stationary  peaks  emergirg 
and  forming  what  appears  to  be  an  almost  isotropic  field.  The  peaks 
probably  are  culturally  related.  Surface- mode  energy  at  NPNT  (lobe 
1)  appears  from  the  general  direction  of  LPC's  1  and  2  and  the  Hudson 
Bay  area.  TFO  has  a  surface-mode  peak  from  the  southwest  (lobe  E), 
generated  either  locally  or  from  the  Baja  coast.  Weak  surface-mode 
energy  persists  from  LPC  14  in  the  northeast  (lobe  14). 


Figure  III-9,  16  October  1964,  F  =  1.  0  Hz 


Regional  storm  energy  again  dominates  the  spectra  at 
1.0  Hz,  generating  surface-mode  energy  at  CPO,  OONW,  NPNT,  and 
GGGR  and  bodywave  energy  at  TFO.  Coherent  power  levels  are  gen¬ 
erally  higher  than  at  0.  6  and  0.8  Hz. 

JLPC  1  is  again  a  source  of  coherent  energy  at  CPO 
(lobe  1).  Also,  bodywave  energy  and  surface-mode  energy  appear  at 
TFO  (lobe  14)  from  the  complex  LPC  14  off  Newfoundland.  A  lesser  peak 
persists  from  the  southwest  (lobe  E).  Surface-mode  energy  is  seen  at 
NPNT  from  the  east  (lobe  14)  and  southwest  (lobe  3),  corresponding 
to  LPC's  14  and  3,  respectively. 

LPC  12  is  seen  as  surface-mode  energy  (lobe  12)  at 
OONW.  The  seemingly  stationary  peak  from  the  southeast  (lobe  M) 
again  appears.  The  source  of  t.-e  weak  surface-mode  peak  from  the 
east  (lobe  U)  is  not  established.  Three  surface-mode  peaks  appear  at 
GGGR:  lobe  V  appears  from  the  west,  possibly  the  result  of  surf 
action  along  the  coast  of  F ranee;  lobe  13  appears  from  the  north,  prob¬ 
ably  generated  in  the  North  Sea  by  LPC's  12  and  13;  and  lobe  A  appears 
from  the  south,  apparently  from  a  stationary  source. 
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Figure  III- 10,  16  October  1964,  F  =  1. 2  Hz 


Coherent  energy  seen  at  1.  2  Hz  does  not  appear  to  be 
primarily  related  to  regional  storm  activity  during  either  noise  sample. 
The  spectrum  for  TFO  is  omitted,  as  unit  cell  effects  prevent  mean¬ 
ingful  interpretation. 

A  surface-mode  peak  from  the  southwest  (lobe  13)  at 
OONW  apparently  is  related  to  LFC  13  in  the  North  Sea.  Weaker  peaks 
from  the  northeast  (lobe  J)  and  southeast  (lobe  M),  also  seen  at  0.8  Hz, 
appear  as  stationary  local  sources.  The  bodywave  peak  at  GGGR  (lobe 
O)  appears  to  be  unrelated  to  meteorological  data;  its  source  is  unde¬ 
fined.  A  surface-mode  peak  from  the  northeast  (lobe  D)  again  appears 
as  an  apparently  stationary  source. 

The  bodywave  peak  at  NPNT  (lobe  P)  is  from  the  general 
direction  of  LPC's  1  and  2.  Also  seen  is  surface-mode  energy  from 
the  southwest  (lobe  3)  related  to  LPC  3.  Coherent  energy  at  CPO  is 
low,  and  the  dominant  coherent  energy  appears  as  an  apparently  sta¬ 
tionary  surface-mode  peak  from  the  southeast  (lobe  B).  Lesser  peaks 
generally  from  the  east  correspond  to  low-level  local  or  coastal  sources, 
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Figure  III- 11,  16  October  1964,  F  =  1. 4  Hz 


A  strong,  very  narrowband  peak  is  observed  from  the 

north-northwest  (lobe  Q)  in  both  data  samples  at  CPO.  This  peak, 

7 

which  has  been  documented  previously,  appears  to  represent  a  local 
stationary  source. 

Bodywave  energy  is  seen  at  NPNT  from  the  southeast 
(lobe  14),  related  to  LPC  14  in  the  Newfoundland  region,  with  a  pos¬ 
sible  contribution  from  LPC  1  along  the  eastern  Atlantic  coast. 

The  apparently  stationary  surface-mode  peak  from 
the  southwest  (lobe  C)  is  the  major  contribution  at  GGGR.  Surface¬ 
mode  energy  also  dominates  at  OONW  (lobe  13),  probably  from  LPC  13. 
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Figure  III- 12,  29  January  1965,  F  =  0.  2  Hz 


Storm -gene rated  energy  again  is  dominant  at  0.  2  Hz 
on  29  January  1965. 

LPC  20  over  Newfoundland  appears  to  be  generating 
surface-mode  energy  at  CPO  and  NPNT  (lobe  20).  TFO,  however, 
shows  only  a  peak  from  the  west  (lobe  26)  at  an  intermediate  velocity. 
This  energy  may  be  related  to  LPC  26  or  may  be  generated  locally. 

The  intense  LPC  22  over  the  Kara  Sea  is  seen  at  OONW 
as  surface-mode  energy  from  the  northeast  (lobe  22).  The  broad, 
poorly  defined  peak  at  GGGR  may  result  either  from  LPC  22  or  from 
the  weaker  LPC  24,  also  located  north  of  Germany. 
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Figure  III- 13,  29  January  1965,  F  =  0,  4  Hz 


The  narrowband  character  of  the  coherent  energy  is 
well-illustrated  by  this  0.  4-Hz  sample.  Ail  stations  report  coherent 
energy  sources  different  from  those  observed  at  0.  2  Hz.  One  possible 
explanation  is  frequency-dependent  earth  filtering,  since  the  stronger 
sources  such  as  LPC's  20  and  22  are  at  considerable  ranges. 

At  CPO,  energy  is  seen  from  the  northwest  (lobe  17), 
apparently  generated  by  the  frontal  passage  and  by  LPC's  16,  17,  and 
18  located  over  the  Plains  states.  Two  peaks  are  suggested,  indica¬ 
ting  both  surface-mode  and  intermediate-velocity  energy.  TFO,  how¬ 
ever,  shows  bodywave  energy  from  the  north  (lobe  19),  apparently 
related  to  LPC  19.  The  peak  at  NPNT  (lobe  22)  has  rather  high  velocity 
for  surface-mode  propagation  but  is  probably  due,  at  least  indirectly, 
to  LPC  22. 

OONW  shows  surface-mode  energy  from  the  northv/est 
(lobe  23),  primarily  from  LPC  23  in  the  Norwegian  Sea.  GGGR  ex¬ 
hibits  the  strong,  narrowband,  apparently  stationary  peak  from  the 
south  (also  reported  on  16  October  1964).  A  weaker  peak  from  the 
north  is  indicated  (lobe  22),  probably  related  to  either  LPC  22  or  24. 
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Figure  III- 14.  High -Re solution  f-k  Spectra 
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Figure  III- 14,  29  January  1965,  F  =  0.  6  Hz 

Again,  at  0.6  Hz  (and  also  at  0.8  Hz),  regional  storm 
activity  ceases  to  be  a  dominant  contributor,  as  seemingly  stationary 
peaks  emerge  at  several  stations. 

A  bodywave  peak  from  the  southwest  (lobe  E)  at  TFO 
corresponds  to  the  general  direction  for  energy  from  LPC  26,  but 
the  indicated  velocity  (15  km/sec)  suggests  a  more  distant  source. 
Bodywave  energy  also  is  seen  at  this  frequency  from  the  southwest  on 
16  October  1964.  This  peak  apparently  is  not  related  to  meteorological 
sources .  Also  indicated  in  both  samples  is  a  weaker  peak  from  the  southeast 
(lobe  F).  A  surface-mode  peak  from  the  southwest  at  CPO  (lobe  B) 
appears  stationary. 

NPNT  shows  a  strong  surface-mode  source  in  the  south 
(lobe  17),  generally  corresponding  to  the  line  of  LPC's  over  the  Great 
Plains.  A  weaker  surface-mode  peak  from  the  north  (lobe  22)  persists 
from  0.4  Hz  and  is  probably  related  to  LPC  22. 

A  surface-mode  peak  from  the  southwest  (lobe  C)  also 
is  seen  on  16  October  1964  at  GGGR  and  is  apparently  from  a  station¬ 
ary  source.  Surface-mode  energy  at  OONVr  from  the  west  (lobe  20) 
may  be  created  either  by  the  intense  LPC  20  over  Newfoundland  or  by 
the  nearer  LPC  23  off  the  Norwegian  coast. 
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Figure  III- 15,  29  January  1965,  F  =  0.  8  Hz 


A  relative  null  in  coherent  energy  level  is  again  generally 
evident  at  0.8  Hz.  Bodywave  energy  appears  at  TFO,  GGGR,  and  OONW; 
surface-mode  energy  is  seen  at  all  stations. 

The  bodywave  peak  from  the  northwest  (lobe  20)  at  OONW 
probably  is  generated  by  LPC  20,  located  over  Newfoundland.  Smaller 
peaks  from  the  northeast  (lobe  J)  and  southeast  (lobe  M)  occur  also  on 
16  October  1964  and  are  apparently  stationary.  Bodywave  energy  is  seen 
also  at  GGGR  from  the  northeast  (lobe  22),  corresponding  to  LPC  22.  A 
surface-mode  peak  from  the  southwest  (lobe  25)  is  seen,  representing 
energy  from  LPC  25  or  the  stationary  source  C. 

The  CPO  spectrum  is  a  virtual  repeat  of  the  one  shown  at 
thie  frequency  for  16  October  1964,  with  seemingly  stationary  peaks  from 
the  northeast  (lobe  K),  southeast  (lobe  B),  southwest  (lobe  G),  and  north 
(lobe  L).  All  are  probably  locally  generated  or  related  to  coastal  surf 
action.  Surface-mode  energy  is  seen  from  the  east  at  NPNT  (lobe  20), 
related  to  LPC  20. 

A  weak  bodywave  peak  is  seen  at  TFO  (lobe  20),  apparently 
from  LPC  20.  Lobe  E  is  generated  by  the  apparently  stationary  surface- 
wave  source  southwest  of  TFO.  Surface-mode  peaks  from  the  north  (lobe 
17)  and  west  (lobe  26)  are  apparently  related  to  LPC's  17  and  26,  respec¬ 
tively.  Surface-mode  energy  also  appears  from  the  northwest,  apparently 
generated  along  the  Pacific  coast  and  related  to  LPC  19.  Due  to  seismom¬ 
eter  spacings  within  the  TFO  crossarray,  the  unit  cell  does  not  extend  to 
this  point  in  the  wavenumber  plane  at  this  frequency;  therefore,  the  peak 
ha 8  space- aliased  to  a  point  within  the  unit  cell  (lobe  19A),  creating  a 
pseudo-source  from  the  southwest. 
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Figure  III- 16,  29  January  1965,  F  =  1 . 0  Hz 


Unlike  the  16  October  1964  sample,  regional  storm 
activity  does  not  dominate  spectra  at  1 . 0  Hz  for  this  sample.  Gen¬ 
erally,  lower  coherent  energy  levels  expose  minor  structural  peaks, 
and  interpretation  is  correspondingly  more  complex. 

For  example,  the  CPO  spectrum  resembles  that  seen 
at  0.8  Hz,  with  the  appearance  of  several  apparently  stationary  sur¬ 
face-mode  peaks.  An  additional  bodywave  peak  from  the  southwest, 
with  unknown  source,  is  seen  (lobe  T).  Two  relatively  high-velocity 
peaks  at  TFO  (lobe  19A)  are  apparently  space-aliased  into  the  unit 
cell  from  the  north  and  correspond  to  coastally  generated  surface¬ 
mode  energy.  Such  interpretation  is  tenuous  and  is  done  on  the  basis 
of  plausible  energy  sources  and  propagation  velocities. 

A  broad  smear  of  bodywave  energy  at  NPNT,  accom¬ 
panied  by  a  surface-mode  peak  (lobe  22),  is  probably  related  to  Arctic 
Ocean  disturbances  created  by  LPC  22.  Low-velocity  peaks  from  the 
east  (lobe  20)  and  south  (lobe  17)  correspond  to  LPC's  20  and  17,  re¬ 
spectively. 

A  surface-mode  peak  (lobe  24)  is  seen  at  GGGR  from 
either  LPC  24  or  the  distant  but  more  intense  LPC  22.  A  much  weaker 
peak  from  the  southeast  (lobe  25)  is  possibly  generated  by  surf  action  from 
LPC  25.  The  indicated  peak  at  240  (2.  0  km/sec)  is  apparently  energy- 

aliased  into  the  unit  cell  from  lobe  24.  Su’  '■»  ce  -  mode  energy  is  also 
seen  at  OONW  (lobe  U)  and  is  possibly  cau6o  by  LPC  24. 
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figure  III- 17,  29  January  1965,  F  =  1. 2  Hz 


Bodywave  peaks  appear  at  both  OONW  and  GGGR  (lobe 
20),  apparently  generated  by  the  intense  LPC  20  over  Newfoundland. 

Surprisingly,  the  stations  located  on  the  same  continent 
as  the  storm  center  fail  to  detect  it.  Instead,  surface-mode  energy 
from  the  southwest  at  NPNT  (lobe  P)  appears  to  be  locally  generated. 

similarly,  an  apparently  stationary  peak  (lobe  B)  south¬ 
west  of  CPO  appears  as  the  only  coherent  structure. 
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Figure  III- 18,  29  January  1965,  F  =  1.  4  Hz 


The  apparently  stationary  narrowband  peak  (lobe  Q) 
is  again  evident  at  CPO  for  1. 4  Hz. 

Two  surface-mode  peaks  appear  at  NPNT:  lobe  22 
from  the  southwest  (probably  related  to  storm  activity  in  the  Arctic 
Ocean)  and  lobe  P  from  the  southeast  (corresponding  to  an  apparently 
stationary  local  source  which  is  possibly  surf  action). 

A  surface-mode  peak  at  OONW  (lobe  24)  is  apparently 
related  to  LPC  24.  Lesser  peaks  from  the  south  (lobe  H)  and  south¬ 
west  (lobe  J)  are  probably  locally  generated,  although  lobe  J  is  pos¬ 
sibly  aliased  into  the  unit  cell  from  the  northeast. 

The  source  of  bodywave  energy  at  GGGR  from  the 
east  (lobe  O)  is  undetermined.  A  similar  peak  is  noted  at  GGGR  for 
1.2  Hz  in  the  16  October  sample. 

Difficulty  is  experienced  in  determining  the  actual  di¬ 
rection  of  the  lesser  surface-mode  peak  apparently  located  on  the 
northwest  edge  of  the  unit  cell. 
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Figure  III-19,  16  October  1964  and  29  January  1965,  F  =  1.6,  1.8,  2.0  Hz 

Spectra  shown  in  this  figure  are  from  CPO  for  both  noiee 
samples  and  for  the  frequencies  1.6,  1.8,  and  2.  0  Hz.  Because  of  space¬ 
aliasing  effects,  CPO  is  the  only  network  station  yielding  interpretable 
f-k  spectra  at  these  frequencies. 

Several  low-level  surface-mode  peaks  appear  at  each  fre¬ 
quency,  generally  representing  a  variety  of  local  sources  that  closely  re¬ 
semble  an  isotropic  field. 

The  only  apparently  storm- related  peaks  are  lobe  1  (16 
October  1964)  and  lobe  17  (29  January  1965). 

The  source  of  the  bodywave  peak  at  1.6  Hz  on  16  October 
1964  is  not  established. 
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SECTION  IV 


BODYWAVE  NOISE  STUDIES 

Several  approaches  to  the  problem  of  partitioning  the  noise 
have  been  examined  during  this  study.  One  approach,  which  uses  high- 
resolution  f-k  spectra  to  examine  the  spatial  distribution  of  the  noise,  is 
discussed  in  Section  III.  Results  of  ''wo  other  approachet  —  a  synthesis 
method  involving  solution  of  simultaneous  equations  and  the  computation 
of  1 -dimensional  wavenumber  spectra  (K-line)  —  are  discussed  in  this 
section. 


The  first  technique  provides  an  estimate  of  the  relative  per¬ 
centages  of  high-  and  low- velocity  energy  for  a  given  noise  sample.  This 
estimate  is  obtained  through  an  analysis  of  spectra  computed  for  the  average 
and  straight-sum  outputs,  along  with  the  array  straight-sum  response.  Cer¬ 
tain  assumptions  are  made  regarding  noise  isotropy  and  random-noise  con¬ 
tent.  A  set  of  simultaneous  equations  is  then  solved  for  the  relative  amounts 
of  high-  and  low-velocity  noise  as  a  function  of  frequency. 

Although  the  results  of  such  a  technique  depend  heavily  on 
the  quality  of  the  assumptions  regarding  the  noise  field  at  a  particular  sta¬ 
tion,  the  assumptions  represent  educated  guesses  from  available  knowledge 
for  that  station.  Undoubtedly,  refinements  to  the  rather  elementary  ap¬ 
proach  would  yield  more  exact  results;  however,  general  conclusions  re¬ 
garding  the  extent  to  which  a  particular  station  is  bodywave -noise-limited 
are  quite  readily  made. 

I 

I 

CPO  and  TFO  are  the  two  stations  investigated  in  this  man¬ 
ner.  Results  indicate  that  CPO  is  primarily  dominated  by  surface-mode 
energy  while  TFO  is  bodywave-noise-limited  except  at  the  very  low  fre¬ 
quencies  (below  0.3  Hz).  These  rather  expected  results  serve  primarily 
to  verify  the  procedure  itself;  however,  the  good  agreement  with  previous 
work  is  a  qualitative  check  on  the  assumptions  made  for  each  station. 
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The  second  technique  involves  analysis  of  K-line  spectra 
computed  for  the  cro.sarray  arm,  a.  the  LRSM  station..  The  CPO  and 
TFO  geometries  are  unsuitable  for  K-line  analysis.  Resolution  is  some- 

monte, ers.  However,  analysis  of  the  spectra  indicate,  predominantly 
sur  ace-mod.  propagation  generally  correlated  with  regional  storm  source. 
None  of  the  three  LRSM  array.  (OONW,  NPNT.  and  GGGR,  appear  to  be 
bodywave- noise -limited,  although  high-velocity  energy  is  a  greater  propor¬ 
tion  of  the  total  a,  NPNT  than  a,  the  other  two  stations. 

likelih  d  a  °th'r  aPPrMCh'S  *°  n°i<e  Partitioning,  such  as  maximum- 
un .  t  h,  ",  m‘er*  a"d  Wiener  dUk-m°del  MCF'*’  «»«  been  found 

«  e  for  vartou.  reasons.  The  adaptive  approach  tend,  to  cause  re- 

1  ction  of  both  bodywave  and  surface-mode  energy.  At  their  present  state 
o  development,  adaptive  di.h-model  technique,  are  not  economically  fea.i- 
bu  for  this  application.  LRSM  array  geometries  generally  have  been  in¬ 
adequate  for  effective  Wiener  disk-model  design. 

An  obvious  requirement  for  meaningful  future  analysis  of 
network  bodywave  noise  is  ,h.  availability  of  high-quality  seismic-array 
data  Since  the  low-v.l„c,,y  energy  doe.  not  appear  coherent  a,  ,h.  net- 
evel,  fUrther  WOrk  ,hould  concentrate  on  studying  the  high-velocity 

componen,.  statlon  array  geometries  ,Wd  be  adequate  proWde  >ig. 

n  Lean,  reduction  in  the  iow-velocity  component  with  station  preprocessing. 

.  developmental  work  continues,  adaptive-filtering  procedure,  will  con¬ 
ceivably  become  suitable  for  bodywave  noise  extraction  a,  the  station  level. 

A.  MULTICHANNEL-FILTER  AND  STRAIGHT-SUM  TECHNIQUES 

1  |fcr  ,  i”  *he  investigation  of  the  relative  quantities  of  high-  and  low- 
velocity  noise  power  contained  in  the  CPO  and  TFO  noise  fields.  ,he  conven¬ 
tional  array  processing  technique,  of  multichannel  filtering  and  straight 
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summing  have  been  utilized  in  a  simplified  method  for  partitioning  a  noise 
field  into  high-  and  low- velocity  energy  components. 

Used  in  the  analysis  are  noise  samples  recorded  at  CPO  on 
16  October  1964  and  29  January  1965  and  at  TFO  on  29  January  1965.  Each 
noiee  sample  is  analyzed,  using  a  6.67-min  time  gate.  Results  of  the  straight 
sum  technique  are  presented  for  each  noise  sample;  however,  results  of  multi¬ 
channel  filtering  are  presented  only  for  CPO  on  16  October  1964. 

A  0. 18-sec  disk  signal-model  MCF  has  been  designed  from 
Whitened  correlations  of  CPO  measured-noise  statistics.  The  disk  signal 
model  has  been  constrained  to  have  an  idealized  gain  of  unity  for  energy  pro¬ 
pagating  at  apparent  horizontal  velocities  of  10  km/sec  or  greater  and  of 
zero  for  lower-velocity  energy.  The  signal  model  has  been  further  constrained 
to  have  tne  shape  of  the  16  October  1964  noise  spectrum  with  an  SNR  of  4. 

Power-density  spectral  estimates  computed  from  ±5.  0-sec 
Parzen-smoothed  correlations  are  presented  for  each  technique  in  Figure 
IV-1.  Also  shown  for  each  noise  sample  are  average  power-density  spec¬ 
tra  formed  from  all  array  inputs.  The  spectral  estimates  are  shown  in  Jb 
relative  to  1  mu  at  1.  0  Hz. 

In  partitioning  the  noise  field,  the  existence  of  only  three 
energy  components  is  assumed:  low  velocity,  P^  high  velocity,  and 
random,  PR  .  P  represents  power.  Low-velocity  energy  is  considered 
to  be  concentrated  in  a  2-km/sec  to  6-km/eec  apparent  horizontal-velocity 
band  and  high-velocity  energy  in  a  8-km/sec,  or  greater,  band.  Random 

energy  is  estimated  to  be  the  total  energy  at  4.  0  Hz  and  is  considered  to  be 
white  across  the  Nyquist  band. 
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For  the  straight  sum,  average  responses  of  K  and  K  are 

L  H 

computed  for  the  low  ■  and  high-velocity  regions,  respectively,  where  K 
represents  that  portion  of  each  component  seen  by  the  straight- sum  response. 
The  energy  in  each  region  is  considered  to  be  uniformly  isotropic.  The  ran¬ 
dom-noise  response  is  considered  to  be  1/N  where  N  is  the  number  of  input 
sensors. 

The  average  input  power  spectrum  may  be  expressed  as  a 
function  of  frequency  by 

PL(f)  +  P„(f)  +  PR  -  PAVG(0  (4-D 

The  straight-sum  spectrum  may  be  expressed  as 

KL(f)  PL(f)  +  KH(f)  PH(f>  +  £  PR  -  PSS(D  (4-2) 


Solution  of  simultaneous  Equations  4-1  and  4-2  for  the  relative 
percentage  of  low-  and  high-velocity  energy  in  the  average  input  spectrum 
yields 


Average  MCF  responses  of  and  are  computed  for  the 
low-  and  high-velocity  regions,  respectively.  Also  computed  is  the  average 
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random-noise  response  of  the  MCF.  The  MCF  power-density  spectrum 
can  then  be  expressed  as 

KL(£)  PL(f>  *  KH(f>  PH(f»  *  KR<°  R  *  PMCF,f>  <4-5> 


Solution  of  simultaneous  Equations  4-1  and  4-4  for  the  relative 


percentage  of  low-  and  high-velocity  energy  in  the  average  input  spectrum 
yields 

(4-6) 


and 


(4-7) 


PH« 


[ 


PMCF(f) 

PAVG(f) 


KL(f,+  P 


AVG 


(f) 


(kl  (f|  -  Kr  <f>j|  '[khIO-KlIO 


Figure  IV-2  shows  the  results  of  applying  each  of  the  two  tech¬ 
niques  to  the  appropriate  noise  samples.  The  illustration  also  presents  a 
composite  straight-sum  result  for  CPO  on  16  October  1964  and  29  January 
1965.  Results  are  evaluated  only  to  1.  5  Hz. 

Results  obtained  for  CPO  on  16  October  1964  for  the  straight- 
sum  and  MCF  techniques  agree  reasonably  well;  both  indicate  .an  average 
low-velocity  component  of  85  percent  and  a  high-velocity  component  of  15 
percent.  Previous  studies  have  indicated  that  the  CPO  noise  field  is  90  per- 

7 

cent  predictable  in  this  frequency  range.  The  noise  field  can  then  be  in¬ 
terpreted  as  consisting  of  76- percent  coherent  low- velocity  energy  and  14- 
percent  coherent  high-velocity  energy.  The  remaining  energy,  then,  is 
due  to  seismic  noise  of  random  nature. 
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Percentage  of  High-  and  Low-Velocity  Noise  Components  at  CPO  and  TFO 


Results  using  the  straight- sum  technique  for  CPO  data  recorded 
on  29  January  1965  indicate  a  significant  minimum  in  the  low-velocity  per¬ 
centage  at  0.  5  Hz.  The  composite  result  indicates  good  agreement  between 
samples,  except  in  the  0.  3-  to  0.  7-Hz  frequency  band.  The  increase  in  low- 
velocity  energy  during  16  October  1964  is  believed  to  be  due  to  intense  storm 
activity  east  of  the  station  on  16  October  (tropical  storm  Isabel)  but  could  be 
due  to  smearing  of  the  0.25-Hz  microseismic  peak  through  smoothing.  If 
such  is  the  case  for  the  differences  observed  between  0.  3  and  0.7  Hz,  the 
overall  ratio  of  high-velocity  to  low-velocity  energy  appears  relatively  con¬ 
stant  above  0.  3  Hz  —  even  in  the  presence  of  increased  microseismic  activ¬ 
ity. 

Results  using  the  straight-sum  technique  for  TFO  data  re¬ 
corded  on  29  January  1965  indicate  that  the  noise  field  is  dominated  by  high- 
velocity  energy  above  0.3  Hz  and  by  low-velocity  energy  below.  This  agrees 
with  previous  studies  of  the  TFO  noise  field  evaluated  by  more  sophisticated 
methods.  ® 

The  2-dimensional  high- resolution  frequency- wavenumber 
(f-k)  spectra  (Section  III)  indicate  no  high-velocity  energy  at  CPO  during  the 
16  October  1964  or  29  January  1965  samples  and  indicate  both  high-  and  low- 
velocitv  energy  at  TFO.  Inherent  in  the  computation  of  these  f-k  spectra  is 
the  assumption  of  planar-wave  propagation.  Possibly  severe  deviations  of 
the  high-velocity  energy  from  planar-wave  propagation  may  exist  due  to 
scattering  and  to  variations  in  trapped- mode  layer  thicknesses,  and  there 
is  the  probability  that  high-velocity  energy  will  be  observed  only  when  it 
constitutes  an  appreciable  portion  of  the  total  energy. 


B.  ANALYSIS  OF  K-LINE  WAVENUMBER  SPECTRA  FROM  LRSM  DATA 


The  following  discusses  an  analysis  of  1 -dimensional  (K-line) 
wavenumber  spectra  obtained  from  6. 67-min  ambient  noise  samples  recorded 
at  LRSM  crossarrays  OONW,  GGGR,  and  NPNT  on  16  October  1964  and  29 
January  1965.  Analysis  objectives  include  a  study  of  noise-power  spatial  dis¬ 
tribution,  comparison  with  high- resolution  2-dimensional  wavenumber  spectra 
computed  for  the  same  samples  (presented  in  Section  III),  and  determination 
of  the  relative  quantities  of  high-  and  low-velocity  noise  power  as  a  function 
of  frequency  at  each  station. 


The  K-line  spectra  are  physically  interpreted  as  projections 
of  2-dimensional  wavenumber  spectra  on  perpendicular  axes  which  are  paral¬ 
lel  to  the  two  crossarray  arms.  Thus,  the  wavenumb  !r  spectra  give  the  power 
density  of  the  ambient  seismic  noise  as  a  function  of  its  apparent  wavenumber 
along  each  of  the  two  arms.  The  basic  input  data  for  calculating  these  wave- 
number  spectra  are  the  cro&spower  matrix  values  ^$.^(f)J  where  f  is  fre¬ 
quency  and  i  and  j  range  over  the  four  seismometers  each  array  arm. 


1.  Description  of  K-Line  Spectra 


To  facilitate  presentation,  K-  line  spectra  for  the  two  noise 
samples  are  shown  (for  each  station)  at  0.  1-Hz  increments  from  0.2  to  1. 0 
Hz,  excluding  0.3  Hz  which,  in  all  cases,  is  very  similar  to  the  0.2-Hz  plot. 

For  each  frequency,  a  wavenumber  power-density  spectrum, 
superimposed  with  an  integrated  wavenumber  power-density  function,  is  pre¬ 
sented  for  each  arm.  For  each  wavenumber  spectrum,  the  0-db  level  indi¬ 
cates  the  spectrum’s  average  value.  The  solid  vertical  line  at  K  *  0  repre¬ 
sents  infinite  apparent  velocity;  i.  e.  ,  the  wave  propagation  is  perpendicular  to 
the  direction  in  which  the  wavenumber  spectrum  is  measured.  Dashed  lines 
on  either  side  of  this  line  show  apparent  velocities  of  8,  3,  and  2  km/sec. 
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The  2-km/sec  line  falls  at  the  extreme  edge  of  the  1.  O-Hz  plots.  Wavenumber 
aliasing  occurs  at  0.  5  cycles/km  for  1-km  spacing  between  seismometers. 

(See  the  appendix  for  array  configurations.)  The  right  halves  of  the  wave- 
number  spectra  show  the  energy  propagating  generally  from  north  to  south, 
while  the  left  halves  represent  energy  moving  from  south  to  north. 


is 


The  normalized  integrated  wavenumber  power-density  function 


r  K 

J  P(k)  dk,  with  f  P(k)  dk  =  1 
K  -K 


This  function  is  shown  plotted  in  fractional  power  vs  wave- 
number  in  cycles/km.  These  integrated  spectral  plots  are  particularly  use¬ 
ful  in  that  they  facilitate  measurement  of  the  amount  of  power  within  any 
given  apparent  wavenumber  (or  velocity)  band. 

2.  Presentation  and  Analysis  of  K-Line  Wavenumber  Spectra 

The  wavenumber  spectra  are  shown  grouped  by  station.  Facing 
each  noise-sample  plot  is  the  related  discussion. 

Some  features  are  common  to  all  plots,  however.  One  such 
feature  is  the  pronounced  drop  in  power  density  outside  the  dashed  lines, 
especially  at  frequencies  below  0. 8  Hz.  This  indicates  the  seismic  validity 
of  the  recorded  data,  since  energy  traveling  at  velocities  lower  than  that  of 
the  fundamental  Rayleigh  wave  (which  is  the  slowest  mode  of  seismic  propa¬ 
gation)  is  not  seismically  interpretable.  Above  0.  8  Hz,  an  appreciable  por¬ 
tion  of  the  power  density  is  apparently  either  random  or  nonseismic;  this  is 
not  particularly  surprising,  since  the  dynamic  range  of  the  FM-recorded 
data  is  virtually  exhausted  at  the  lower  power  densities  observed  for  these 
frequencies. 
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The  velocity  of  the  Rayleigh-wave  energy  appears  to  vary 
with  station  and  ranges  from  about  3.  1  km/sec  at  OONW  to  about  2.  7  km/ 
sec  and  2.4  km/sec  at  NPNT  and  GGGR,  respectively .  Some  dependence 
of  propagation  velocity  with  azimuth  is  noted  at  OONW,  suggesting  a  local 

source.  Also,  there  are  indications  of  surface- wave  dispersion  at  OONW 
and  GGGR. 

Another  striking  feature  is  the  dominance  of  directional 
Rayleigh-wave  energy  at  all  stations,  especially  at  the  lower  frequencies. 
Above  0.  4  Hz,  there  is  some  evidence  of  mantle  P-wave  energy  at  NPNT. 

The  organized  noise  appears  generally  correlated  with  regional  storm  ac¬ 
tivity. 

General  agreement  in  spatial  organization  is  seen  between 
the  K-line  and  the  2-dimensional  spectra  (presented  in  Section  HI)S  especially 
at  the  peak  power  frequency.  Observed  at  other  frequencies  are  some  dif¬ 
ferences  caused  primarily  be  two  effects.  The  first  is  a  smearing  of  spec¬ 
tral  power  across  frequencies,  which  occurs  when  the  K-line  autospectra 
are  boxcar- smoothed  to  stabilize  the  spectral  estimate;  the  2- dimensional 
spectra  are  unsmoothed,  incurring  a  somewhat  higher  variance  of  spectral 
estimate  but  x’etaining  an  extremely  narrowband  look  at  the  noise  field.  The 
second  effect  is  a  general  lack  of  spatial  resolution  in  the  K-line  plots,  which 
is  due  primarily  to  the  availability  of  only  four  seismometers  along  each  arm. 
Burg  has  shown  that  interpretation  using  only  four  seismometers  is  practical 

and  has  described  the  effect  on  spectral  resolution  as  a  function  of  the  number 

8 

•  of  sensors  along  the  arm. 
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Figure  IV-3.  Wavenumber  and  Integrated  Spectra  from  Two 
Arms  at  C'ONW,  29  January  1965 
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Figure  IV-3,  OONW,  29  January  1965 


Directional  surface-mode  energy  dominates  the  K-line  spectra 
during  this  sample. 

At  the  point  of  peak  power,  represented  by  0.  2  Hz  (and  0.  3 
Hz  which  is  not  shown),  there  is  excellent  agreement  with  the  2-dimensional 
spectra  Both  show  Rayleigh- wave  energy  (3.  3  to  3.  5  km/sec),  apparently 
generated  by  LPC  22  over  the  Kara  Sea  northeast  of  OONW.  At  0.4  and  0.  5 
Hz,  a  shift  in  source  direction  is  noted,  as  the  K-line  spectra  indicate  low- 
velocity  (3.  l-km/sec)  energy  from  the  southeast,  probably  caused  by  LPC 
24  over  the  Baltic  Sea.  Energy  at  about  5  km/ sec  from  the  northwest  also 
appears  and  may  be  related  to  LPC  23  over  the  Norwegian  Sea  and  to  surf 
action  along  the  northwestern  coast  of  Norway.  This  latter  source  is  veri¬ 
fied  by  the  2-dimensional  spectra  at  0.  4  Hz,  although  the  observed  velocity 
is  much  lower  (2.  0  km/sec). 

LPC  24  still  appears  as  a  contributor  at  0.6  and  0.7  Hz,  with 
an  added  indication  of  energy  from  the  west,  traveling  at  about  5  and  3  km/ 
sec,  respectively,  at  these  frequencies.  Surface-mode  energy  from  this 
direction  (at  3  km/sec)  is  detected  with  the  2-dimensional  spectra  at  0.6  Hz, 
with  a  bodywave  peak  (8  km/sec)  at  0.8  Hz  from  the  same  azimuth.  On  the 
strength  of  the  bodywave  peak,  LPC  20  (over  Labrador)  is  the  likely  con¬ 
tributor  of  this  energy.  At  0.8,  0.9,  and  1.0  Hz,  the  K-line  spectra  appear 
virtually  isotropic,  and  directional  source  determination  is  not  attempted. 

There  is  indication  of  possible  surface-wave  dispersion.  Re¬ 
sults,  while  not  conclusive,  indicate  reduction  in  horizontal  propagation  ve¬ 
locity  with  increasing  frequency,  ranging  from  about  3.4  km/sec  at  0.2  and 
0.  3  Hz  to  3. 0  or  less  above  0.  5  Hz. 
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Figure  IV-4,  OONW,  16  October  1964 


Surface- mode  energy  again  appears  dominant.  Good  agree¬ 
ment  is  seen  with  the  2-dimensional  spectra  at  0.2  Hz,  with  surface-mode 
energy  (3  km/sec)  from  the  southwest  and  probably  related  to  LPC  13.  At 
0.4  Hz,  the  source  (still  probably  related  to  LPC  13)  appears  to  be  located 
further  south,  with  energy  propagating  at  3.  2  km/ sec.  The  f-k  spectra  in¬ 
dicate  energy  from  this  direction:  however,  seen  in  the  K-line  spectra  at 
8  km/sec  from  the  northwest  is  an  additional  peak  possibly  caused  by  LPC  12. 

At  0.  5,  0.  6,  and  0.  7  Hz,  surface-mode  energy  is  seen  from 
due  south  at  2.4,  2.2,  and  2.7  km/sec,  respectively,  and  agrees  well  with 
the  2-dimensional  spectra  at  0.6  Hz.  The  low  velocity  suggests  a  nearby 
source,  such  as  surf  action  in  the  bay  south  of  Oslo,  which  is  likely  an  in¬ 
direct  result  of  LPC  13.  The  energy  from  LPC  13  appears  to  persist  at 
0.8  and  0.9  Hz,  with  an  additional  high-velocity  source  indicated  from  the 
southeast  (lobe  M)  which  does  not  appear  directly  attributable  to  weather. 

Both  peaks  are  seen  in  the  2-dimensional  spectra  at  0.8  Hz,  although  lobe 
M  appears  as  a  surface-wave  peak.  The  K-line  spectra  again  appear  iso¬ 
tropic  at  1.  0  Hz. 

Slight  evidence  of  dispersion  again  appears,  although  it  is 
somewhat  obscured  by  unusually  low-velocity  energy  from  the  south. 
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Figure  IV-5,  NPNT,  29  January  1965 

Excellent  agreement  between  the  K-line  and  2-dimensional 
spectra  is  seen  at  0.  2  (and  0.  3)  Hz. 

Rayleigh-mode  energy  from  the  southwest  (at  2.7  km/sec) 
apparently  is  related  to  LPC  20  over  Newfoundland.  At  0.4  and  0.  5  Hz,  a 
new  source  is  indicated;  energy  from  the  south- southeast  is  probably  related 
to  the  frontal  passage  over  the  Great  Plains  and  the  associated  LPC'*  16,  17, 
and  18. 

Both  surface-mode  (3.  1  km/sec  at  0.4  and  0.5  Hz)  and  body- 
wave  (8  km/sec  at  0.  5  Hz)  peaks  are  observed.  The  2-dimensional  spectra 
at  0.4  Hz  show  only  a  Rayleigh-mode  peak  from  the  northwest,  which  is  pos¬ 
sibly  related  to  LPC  22.  Additional  energy  is  indicated  from  LPC  20  in  the 
K-line  spectra  at  0.  4  Hz. 

Both  spectra  indicate  surface-mode  energy  from  LPC's  16, 

17,  and/or  18  at  0.6  Hz.  High-velocity  energy  also  is  indicated.  Surface¬ 
mode  energy  from  this  source  is  still  apparent  at  0.7  Hz  in  the  K-line  spec- 
tra.  From  0.8  to  1.  0  Hz,  the  K-line  spectra  do  not  appear  seismically  in¬ 
terpretable. 
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Figure  IV-6,  NPNT,  16  October  1964 

The  K-line  spectra  for  this  sample  indicate  a  rather  dominant 
contribution  from  LPC  14  off  Labrador  (southeast  of  NPNT). 

At  0.2  Hz,  there  is  indication  of  surface-mode  energy  (about 
4  km/  sec)  from  the  southwest,  corresponding  roughly  to  the  2-km/sec  peak 
seen  in  the  2-dimensional  spectra.  Surface-mode  energy  is  observed  from 
LPC  14  at  3. 3,  3.2,  2.8,  and  3. .  km/ sec  in  the  K-line  spectra  for  0.  2,  0.3, 
0.4,  and  0.  5  Hz,  respectively,  and  agrees  with  the  2-dimensionai  peak  seen 
at  0.4  Hz.  However,  bodywave  peaks  are  seen  at  0.6  and  0.7  Hz  in  the  K-line 
spectr"  (also  corresponding  to  LPC  14)  with  no  evidence  of  the  bodywave  peak 

from  the  west  which  dominates  the  2-dimensional  spectra  and  is  apparently 
related  to  LPC  5. 

From  0.  8  to  1.  0  Hz,  bodywave  energy  from  the  southwest, 
po, Bibly  related  either  to  LPC  1  or  LPC  14,  i.  indicated  in  the  K-line  epec- 

tra.  The  2-dimen.ional  epectra  indicate  eurface-mode  energy  from  the  eame 
directions  at  0.  8  and  1.  0  Hz. 

Most  notable  in  this  sample  are  the  predominant  directionality 
of  the  coherent  noise  and  the  large  amounts  of  bodywave  energy,  all  apparently 
storm- generated  and  traversing  the  continent  across  the  Canadian  Shield. 


Figure  IV- 7 ,  GGGR,  29 January  1965 


A  poorly  defined  peak  from  the  north  in  the  2-dimensional 
spectra  at  0.2  Hz  is  seen  on  the  K-line  spectra  as  two  surface-mode  peaks, 
one  traveling  northwest  and  the  other  southeast.  These  probably  are  not  seis- 
mically  valid.  The  rather  peculiar  shape  of  the  GGGR  spectrum  (Section  III) 
should  be  mentioned. 

From  0.3  to  0.7  Hz,  the  K-line  spectra  are  dominated  by 
Rayleigh-mode  energy  from  the  southwest,  which  is  possibly  generated  by 
LPC  25  but  is  more  likely  from  a  stationary  source  since  broadband  energy 
is  also  seen  from  this  direction  on  16  October  1964.  Above  0.6  Hz,  a  per¬ 
sistent  high-velocity  peak  is  also  noted;  whether  this  is  instrumental  or  body- 
wave  energy  from  LPC  20  (northeast  of  GGGR)  is  not  established. 

The  narrowband  stationary  source  in  the  south  (lobe  A),  in¬ 
dicated  in  the  2-dimensional  spectra  at  0.  4  Hz,  is  not  seen  in  the  K-line 
spectra.  Energy  from  LPC  25  is  seen  in  the  2-dimensional  spectra  at  0.6, 
0.8,  and  1. 0  Hz. 

There  is  some  indication  of  surface-mode  dispersion,  with 
the  velocity  of  the  Rayleigh  wave  from  the  southwest  decreasing  from  about 
3.8  km/sec  at  0.3  Hz  (not  shown)  to  1.9  km/sec  at  0.7  Hz.  This  effect  is 
established  also  for  the  16  October  1964  sample. 
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Figure  IV-8.  Wavenumber  and  Integrated  Spectra  from  Two 
Arm*  at  GGGR»  16  October  1964 


Figure  IV-8,  GGGR,  16  October  1964 


GGGR  again  receives  broadband  surface-mode  energy  from 

the  southwest. 

Little  agreement  is  seen  between  the  K-line  and  2-dimensional 
spectra,  except  at  0.  6  Hz. 

At  0.2  and  C.8  Hz,  the  f-k  spectra  indicate  surface-mode 
energy  from  the  northwest  (LPC  13)  and,  at  0.  4  Hz,  from  stationary 
source  A. 

Above  0.8  Hz,  the  K-line  spectra  are  not  seismically  in¬ 
terpretable  (in  terms  of  directional  noise).  From  0.  2  to  0.  7  Hz,  the 
surface-wave  velocity  is  again  observed  o  decrease  as  frequency  increases . 
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C.  ANALYSIS  OF  BODY  WAVE  NOISE  POWER  DENSITY 

Estimates  of  the  bodywave  noise  levels  at  the  five  array  sta¬ 
tions  were  made  for  the  16  October  1964  and  29  January  1965  noise  samples 
(excluding  TFO  for  16  October  1964).  Estimates  for  OONW,  GGGR,  and 
NPNT  were  made  from  joint  analysis  of  the  2-dimen3ional  and  K-line  spectra 
and>  in  particular,  the  integrated  fractional  power  curves  for  each  array  arm. 
Estimates  for  TFO  and  CPO  were  based  on  the  analysis  discussed  in  subsec¬ 
tion  A  of  this  section. 

Figure  IV-9  displays  short-period  absolute  power-density 
2 

spectra  (relative  to  1  m|i  /Hz  ground  motion  at  1  Hz  for  each  station)  for 
the  two  noise  samples.  Bodywave  noise  power- density  levels  also  are  shown 
for  each  station.  Instrument  responses  have  been  removed  to  permit  com¬ 
parison  of  power  levels  across  the  network,  since  different  instruments  were 
in  use  at  the  various  stations. 

Stations  CPO,  OONW ,  and  GGGR  appear  to  be  dominated  by 
low-velocity  energy,  while  TFO  is  apparently  bodywave-noise-limited. 

NPNT  appears  rather  intermediate,  with  the  high-velocity  energy  being  only 
3  to  4  db  below  the  total  power  level  at  frequencies  approaching  1. 0  Hz. 
Bodywave  levels  vary  considerably  from  sample,  to  sample  and  actually  appear 
less  consistent  than  the  total  power  noise  levels.  Such  behavior  likely  re¬ 
flects,  in  part,  a  degree  of  variability  inherent  in  the  analysis  approach  but 
also  suggests  that  the  bodywave  noise  is  probably  not  time- stationary . 

Figure  IV- 10  combines  the  high-velocity  noise  power-density 
plots  from  each  station  for  the  two  noise  samples  to  permit  comparisons  of 
bodywave  noise  levels  across  the  network.  Also  shown  are  the  combined  total 
power-density  plots. 
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Figure  IV- 10.  Comparison  of  Station  Total  and  High-Velocity 
Power-Density  Spectra  for  Each  Noise  Sample 
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Two  characteristics  are  noted:  first,  there  appears  to  be 
considerable  variation  in  levels  among  stations,  suggesting  a  lack  of  space- 
stationarity;  second,  significantly  different  levels  are  seen  for  the  stations 
located  on  the  European  continent  (OONW,  GGGR)  compared  to  the  levels 
for  the  stations  on  the  North  American  continent  (CPO,  TFO,  NPNT).  The 
variation  in  bodywave  spectral  level  between  continents  actually  appears  to 
increase  with  increasing  frequency  up  to  1.0  Hz. 

Because  of  limitations  in  the  procedures  for  estimating  the 
bodywave  noise  components,  the  estimates  are  probably  subject  to  consider¬ 
able  error,  especially  at  the  lower  frequencies  (near  0.  2  Hz).  Errors  ex¬ 
ceeding  3  db  seem  rather  unlikely,  however,  and  will  probably  not  account 
for  the  observed  differences  in  continental  levels.  Several  alternative  ex¬ 
planations  for  this  effect  are  possible:  e.  g.  ,  systematic  error  in  procedure 
including  calibration  of  power  spectra  and  instrument  response  removal  or 
instrumentally  generated  (random)  noise  levels  induced  by  differences  in 
the  dynamic  ranges  of  the  recordings.  More  likely,  however,  the  effect  is 
real  and  is  related  to  differences  in  area  sensitivities  to  bodywave  energy. 
Such  differences  may  arise  from  geological  structure  or  relative  distances 
to  P-wave  sources  such  as  storms  and  seismically  active  areas. 
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APPENDIX 


DESCRIPTIONS  OF  WORLDWIDE  NETWORK 
AND  DATA  LIBRARY 

Figure  A-l  shows  the  relative  locations  of  the  10  seismologies! 
data  recording  stations  forming  the  network  displayed  on  a  Mercator  projection 
of  the  world.  Nine  of  the  stations  are  located  in  the  northern  hemisphere, 
and  one  (LZBV)  is  in  the  southern  hemisphere.  The  stations  were  selected 
to  provide  simultaneous  data  recordings  over  a.  extensive  an  area  as  possible. 
Table  A-l  presents  pertinent  information  regarding  each  station;  i.  e. ,  name, 
location,  elevation,  number  of  short-period  instrumeofs,  and  predominant 
basement  rock. 


Table  A-l 

GENERAL  STATION  INFORMATION 


Station 

Designator 

Station  Location 

j  Station  Coordinates 

Instruments 

Latitude 

n 

Elevation 

(km) 

Type* 

No. 

Composition 

ADIS 

Adak,  Aleutian  Islands 

51.87  N 

176.67  W 

0.  06 

BEN 

1 

Basalt 

CPO 

McMinnville,  Tennesse 

35,  58  N 

85.  57  W 

0.  57 

JM 

10 

Limestone 

DHNY 

Delhi,  New  York 

42.23  K 

74.  88  W 

0*65 

BEN 

1 

Sandstone 

GGGR 

Grafenburg,  West  Germany 

49.  6BN 

11.20  E 

0.  53 

BEN 

I 

Limestone 

HWIS 

Hawaiian  Islands 

19.  97  N 

155.70  W 

0.71 

BEN 

Basalt 

LZBV 

La  Pae,  Bolivia 

16.25  S 

68.47  W 

3.99 

JM 

H 

Limestone 

NPNT 

Mould  Bay,  N,  W.  Territory 

76.25  N 

119.37  W 

0.  06 

JM 

Alluvium 

OONW 

Oslo,  Norway 

61.  05  N 

10.85  E 

0.  56 

BEN 

Glacial  till 

RKON 

Red  Lake,  Ontario 

50.  83  N 

93. 67  W 

0.37 

BEN 

Granite 

TFO 

Payson,  Arieona 

34.27  N 

111.27  W 

1.49 

JM 

Granite 

BEN:  Short-period  vertical  Benioff 

JM  :  Short- period  vertical  Johnaon-Matheaon 

All  atatlona  uae  3-component  Sprengnethe?  for  long-period  recording 


Each  of  the  stations,  except  CPO  and  TFO  which  are  VELA 
observatories,  are  Geotechnical  Corporation  LRSM  mobil  installations.  Seven 
of  the  10  stations  have  short-period  arrays.  Configurations  for  each  of  the 
array  stations  are  shown  in  Figure  A -2. 

Three -component  short-period  and  long-period  instruments 
were  available  at  each  station;  however,  only  vertical  instruments  were  uti¬ 
lized  for  the  bodywave  noise  analysis.  The  short-period  vertical  recordings 
at  each  station  were  taken  from  either  a  Johnson-Mathe son  or  a  large  Benioff 
seismometer  system  (Table  A-l),  while  long-period  recordings  were  made 
with  Sprengnether  instruments.  Responses  for  the  appropriate  seismometer 
systems  are  shown  in  Figure  A-3.  These  response  curves  represent  the 
normalized  magnification  of  a  seismometer  system  output  for  a  constant  cur¬ 
rent  input.  Short-period  curves  are  normalized  at  1.  0  Hz,  long-period  at 
0.  04  Hz. 

Significant  geological  aid  environmental  differences  exist 
between  stations  (Table  A-l),  as  would  be  expected  for  a  world  array.  Geo¬ 
logical  differences  resulf  in  a  variation  of  crustal  filtering  effects,  scattered 
energy  components,  ant  inherent  time  residuals  for  teleseismic  energy.  Varia¬ 
tions  in  the  ambient  noise  levels  at  the  stations  result  from  environmental 
changes  between  stations  such  as  prevailing  winds,  weather  conditions,  surf 
action  along  nearby  coastlines,  and  local  cultural  sources. 

Field  data  at  each  of  the  network  stations  were  recorded  on 
frequency-modulated  (FM)  magnetic  tapes.  The  Seismic  Data  l  aboratory 
(SDL)  subsequently  digitized  and  demultiplexed  these  FM  tapes.  Figure  A-4 
is  a  block  diagram  showing  the  sequence  of  stepB  involved  in  the  data  prepro¬ 
cessing. 
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Figure  A-4.  Flow  Diagram  of  Data  Preprocessing  Sequence 

Data  were  available  at  SDL  in  one  of  two  possible  forms: 
standard  FM  field  tapes  or  compressed  FM  tapes  (which  have  been  rerecorded 
in  comprensed  form  from  standard  field  tapes  to  conserve  data  storage  space). 
The  short-period  data  were  digitised  at  10  samples/ sec  in  a  12 -channel  multi¬ 
plexed  form.  Due  to  the  multiplexing,  an  inherent  sampling  delay  between 
adjacent  channels  of  1/600  sec  and  1/120  sec  was  introduced  from  standard 
and  compressed  FM  tapes,  respectively. 

Prior  to  digitization,  an  analog  bandpass  filter  with  half- 
power  points  at  0.  017  Hz  and  5.  00  Hz  was  applied  to  the  short-period  data. 
The  filter  gain  was  unity  within  ±  1/2  db  inside  the  passband,  and  attenuation 
was  at  24  db/octave  outside  the  passband.  The  long-period  data  were  digi¬ 
tized  at  2  samples/sec  subsequent  to  antialias  filtering  with  a  0.  002-  to  0.25- 
Hz  bandpass  filter.  The  filter  response  roll-off  rates  were  again  approxi¬ 
mately  24  db/octave. 

The  data  after  digitization  were  demultiplexed  by  SDL  into 
IBM -compatible  tapes,  which  were  subsequently  reformatted  at  TI  for  com¬ 
patibility  with  existing  software. 
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Noise  samples  were  selected  for  the  network  on  the  basis  of 
usable  seismic  data  recorded  simultaneously  at  each  station  over  time  inter¬ 
vals  of  approximately  25  min  (short-period)  and  133  min  (long-period).  Data 

$ 

were  chosen  from  35-mm  Develocorder  film,  and  requests  were  submitted 
to  SDL  for  digitization  and  demultiplexing.  Calibration  data  and  station- 
calibration  logs  were  also  obtained  for  the  selected  data.  Table  A-2  describes 
the  noise  ensemble  received. 

Five  noise  samples  were  selected  for  the  short-period  noise 
ensemble.  Two  (16  October  1964  and  29  January  1965)  contained  significant 
storm  activity  and  were  chosen  from  the  ensemble  for  use  in  the  bodywave 
noise  analysis.  Each  was  represented  at  array  stations  CPO,  GGGR,  NPNT, 
OONW,  and  TFO.  All  sensors  were  available  at  each  station  with  the  exception 
of  Z5  at  CPO  and  Z7  at  NPNT  on  16  October  1964. 

Long -period  data  were  obtained  for  the  time  periods  1545  to 
1758  GMT  on  16  October  1964  and  0515  to  0728  GMT  on  29  January  1965.  Sta¬ 
tions  ADIS,  CPO,  RKON,  and  TFO  were  represented  for  both  samples, 

HWIS  for  the  16  October  1964  sample  only,  arid  NPNT  for  the  29  January  1965 
sample  only.  Data  also  were  received  for  OONW  and  NPNT  for  16  October 
1964;  however,  difficulty  was  encountered  in  reformatting  the  SDL  tape  records. 
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Table  A -2 


SHORT -PERIOD  NOISE  ENSEMBLE 


* 

Noise  Sample 

Station 

12  Aug  1964 
1544  to  1610 
(GMT) 

16  Oct  1964 
1630  to  1656 
(GMT) 

6  Nov  1964 
2302  to  2328 
(GMT) 

23  Nov  1964 
0600  to  0626 
(GMT) 

29  Jan  1965 
0600  to  0626 
(GMT) 

ADIS 

+ 

+ 

+ 

+ 

+ 

CPO 

- 

+ 

+ 

+ 

DHNY 

+ 

+ 

+ 

+ 

+ 

GGGR 

+ 

+ 

+ 

- 

+ 

HWIS 

- 

- 

I 

I 

I 

LZBV 

+ 

- 

+ 

+ 

+ 

NPNT 

+ 

+ 

+ 

+ 

+ 

OONW 

+ 

+ 

+ 

+ 

RKON 

+ 

f 

+ 

+ 

+ 

TFO 

- 

+ 

+ 

+ 

+ 
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